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Abstract. We consider a periodic Schrodinger operator and the composite Wannier func- 
tions corresponding to a relevant family of its Bloch bands, separated by a gap from the rest 
of the spectrum. We study the associated localization functional introduced in |MaVaj and 
we prove some results about the existence and exponential localization of its minimizers, in 
dimension d < 3. The proof exploits ideas and methods from the theory of harmonic maps 
between Riemannian manifolds. 
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1. Introduction 

Many transport properties of electrons in crystalline solids are understood by the analysis 
of Schrodinger operators in the form 

(1.1) H = -A + Vr(x) acting in L 2 (R d ) 

where the function Vr : R d — > R is periodic with respect to a Bravais lattice T ~ Z, d . The func- 
tion Vr represents (in Rydberg units) the electrostatic potential experienced by a test electron 
and generated by the ionic cores of the crystalline solid and, in a mean-field approximation, 
by the remaining electrons. We refer to |CaLeLi| ICaDeLe] for the mathematical justification 
of such a model in the reduced Hartree-Fock approximation. 

A crucial problem in solid state physics is the construction of an orthonormal basis, canoni- 
cally associated to the operator H, consisting of functions which are exponentially localized in 
space. Indeed, such a basis allows to develop computational methods which scale linearly with 
the system size |Go| . makes possible the description of the dynamics by tight-binding effective 
Hamiltonians |NaMa[ IHS| . and plays a prominent role in the modern theories of macroscopic 
polarization |KSV1 [Re] and of orbital magnetization in crystalline solids |TCVR| . 

A convenient basis has been proposed by Wannier |Wa] . and Wannier functions (see Defi- 
nition [T]) are nowadays a fundamental tool in solid state physics. The problem of proving the 
existence of exponentially localized Wannier functions was raised in 1959 by the Nobel Prize 
winner W. Kohn |Koj . who solved it in dimension d = 1 in the case of a single isolated Bloch 
band for a centrosymmetric potential. The latter condition has been later removed by J. des 
Cloizeaux |C12| . In higher dimension, the problem has been solved, always in the case of a 
single isolated Bloch band, by J. des Cloizeaux |C11| IC12] for centrosymmetric potentials and 
finally by G. Nenciu under general hypotheses (Nelj, see also |HS| for a different proof and 
the review paper |Ne2| . As for dimension d = 1, an alternative approach based on the band 
position operator has been developed, yielding exponential localization of Wannier functions 
even for non periodic Schrodinger operators, see |NN[ ICNN| and references therein. 

However, in dimension d > 1 the Bloch bands of crystalline solids are not, in general, 
isolated. Thus the interesting problem, in view of real applications, concerns the case of 
multiband systems, and in this context the more general notion of composite Wannier functions 
is relevant |B1[ IC11] . The existence of exponentially localized composite Wannier functions has 
been proved in |Ne2| in dimension d = 1. As for d > 1, this problem remained unsolved until 
recently |Pa[ IBPCM| . when an existence result was obtained by geometric methods. 

To circumvent such conceptual difficulty, and in view of the application to numerical sim- 
ulations, the solid-state physics community preferred to introduce the alternative notion of 
maximally localized Wannier functions [MaVaj . The latter are defined as the minimizers of a 
suitable localization functional, known as the Marzari-Vanderbilt (MV) functional. In |MaVa| 
it is also conjectured that the minimizers, whenever they exist, are exponentially localized. 
While such an approach provided excellent results from the numerical viewpoint, being nowa- 
days an every-day tool in computational physics |MYSV| IMSVj . a mathematical analysis of 
the MV functional is still missing. Our goal is to fill this gap. 

In this paper, we prove preliminarily that minimizers of the MV functional do exist for 
d < 3 in a suitable function space. More relevantly, we prove that the minimizers of the MV 
functional are exponentially localized. More precisely, let m > 1 be the number of Bloch 
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bands corresponding to the composite Wannier functions; then the result is proved in three 
cases: if m = 1 for any d > 1, if 1 < d < 2 for any m > 1, and if d = 3 under the constraint 
2 < m < 3 (Theorem 16. ip . The proofs are dimension-dependent. In the first two cases, 
exponential localization holds true for any stationary point of the MV functional, i.e. for any 
solution of the corresponding Euler Lagrange equations, while in the three dimensional case 
the minimality property is crucial. In view of that, when d = 3 the constraint on m arises from 
the methods we exploit in the proof; we believe that it is purely technical, since we do not see 
any physical reason which prevents the result from being true for any m > 1. So far, even if 
one is merely interested in almost-exponential localization, i.e. to prove that the maximally 
localized composite Wannier functions decrease faster than the inverse of any polynomial, the 
constraint on m is needed in our approach. 

The existence of a basis of Wannier functions which are exponentially localized in space 
is also relevant for the periodic Pauli and Dirac operators. The Z d -symmetry of the latter 
operators allows a direct integral decomposition analogous to the Bloch-Floquet transform, 
yielding a structure very similar to the one appearing for periodic Schrodinger operators. With 
few modifications, our methods can also be used to discuss the corresponding minimization 
problem in this context. 

Mathematically, the MV functional can be rewritten, after Bloch-Floquet transform and by 
exploiting the time-reversal symmetry of the operator (|l.ip (i.e. the fact that H commutes 
with the complex-conjugation operator), as a perturbation Fmv of the Dirichlet energy for 
maps from T* to U(m), see (TO) , where T* ~ M d /(27rZ) d is a d-dimensional flat torus and 
U(m) C M m (C) is the unitary group. The exponential localization of a minimizer of the MV 
functional is related to the analyticity of the corresponding minimizer of Fmv- The existence 
of a minimizer of the latter functional follows essentially from the direct method of calculus 
of variations. 

We prove the analyticity of the minimizers of -Fmv by adapting ideas and methods from the 
regularity theory for harmonic maps, see |CWY[ ILW1[ IB!] and references therein. The crucial 
step is to prove that any minimizer of Fmv is continuous. In the two dimensional case, this fact 
is a consequence of the hidden structure of the nonlinear terms in the Euler Lagrange equation 
for the Fmv functional. In the three dimensional case, the continuity follows instead from the 
deeper fact that minimizers at smaller and smaller scales look like minimizing harmonic maps 
from to U{m). We are able to prove that, for m < 3, the latter are actually real-analytic 
(Theorem IA.11|) . by showing constancy of the tangent maps as in the important paper |SU2| . 
As a consequence, we obtain the continuity of the minimizers of Fviv and, in turn, analytic 
regularity. 

We hope that our result will contribute to a fruitful exchange of ideas, problems and methods 
between the solid-state physics and the mathematics community, in the study of Schrodinger 
operators and related physical phenomena. 

Acknowledgments. We are grateful to Ch. Brouder, G. Dell'Antonio and S. Teufel for a 
careful reading of the manuscript and for useful comments. We are indebted to A. D'Andrea, 
D. Fiorenza, G. Mondello and P. Papi for stimulating discussions about the geometric and 
Lie-theoretic aspects of the problem. 
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2. WANNIER FUNCTIONS AND Bloch BUNDLES 

2.1. Periodic Schrodinger operators and Bloch-Floquet transform. In this section 
we recall, following [PST] . that, in view of their invariance under T-translations, Schrodinger 
operators in the form (jl.ip can be decomposed as a direct integral of simpler operators by the 
(modified) Bloch-Floquet transform. A comparison with the formalism of the classical Bloch- 
Floquet transform, appearing in most of the physics literature, is summarized in Remark l2.6l 

The lattice T, corresponding to the Bravais lattice in physics, is described as 

r = 1 7 G R d : 7 = X^^=i n i 7j f° r some % G z| , 

where {7i,---,7d} are fixed linearly independent vectors in R d . The dual lattice, with re- 
spect to the ordinary inner product, is V* := {k G R d : k • 7 G 2vrZ for all 7 G T} . To fix the 
notation, we denote by Y the centered fundamental domain of V, namely 

Y = {x£R d :x = Zj=l a 3 H for a j G hi 5]} ■ 
Analogously, we define the centered fundamental domain Y* of T* by setting 
Y* = [k G R d : fc = E -=i^ 7* for fcj- G \—\, 5]} , 

where {7|} is the dual basis to {7^}, i.e. 7* • 7i = 2n5j t i. When the opposite faces of Y* 
are identified, one obtains the torus := M^/r*, equipped with the flat Riemannian metric 
induced by R d . 

For ip G 5(M a! ), one defines the modified Bloch-Floquet transform as 
(2.1) (U BF ip)(k,y) := 5>- ifc '^V(y + 7), y G R d , k G R d . 

One immediately reads the periodicity properties 

(^bfVO (fe, 2/ + 7) = (^bfVO (k, y) for all 7 G Y , 

(2.2) 

(WBF^)(fc + A,y) = e~ iX -y (U BF iP)(k,y) for all A G T* . 

For any fixed k G ffi d , (Z^bfVO ') * s a F-periodic function and can thus be regarded as an 
element of T~L{ := L 2 (T d ), T d being the flat torus M rf /T. On the other hand, the second equation 
in (12. 2ft can be read as a pseudoperiodicity property, involving a unitary representation of the 
group r*, given by 

(2.3) t:T* ^U(U { ), A^r(A), (r(A)p) (y) = e iA ' 
Following |PST| . it is convenient to introduce the Hilbert space 

(2.4) U T ■= G L? oc (M d ,ftf) : tp(k - A) = r(A)p(fc) VA G T* , for a.e. fe G , 
equipped with the inner product 

^)h t = / dk(ip(k), ip(k)) Hf . 
Jy* 



BLOCH BUNDLES, MARZARI-VANDERBILT FUNCTIONAL AND WANNIER FUNCTIONS 5 

Obviously, there is a natural isomorphism between 7i T and L 2 (Y* given by restriction 
from R d to Y*. The map defined by (|2.ip extends to a unitary operator 

Ubf ■ L 2 {R d ) — > U T ~ / H f dlfe, 

Jy* 

with inverse given by 

(U^ip) (x) = [ dk e ifc [x]), 

V ' \Y*\2 JY* 

where [•] refers to the a.e. unique decomposition x = j x + [x], with j x € T and [x] G 1". 
Finally, from the definition (|2. 1 j) one easily checks that 

?/> G W m ' 2 {R d ), m G N ZY BF V> G £ 2 (^*, VF m ' 2 (T d )), 

(2.5) 

(x) m V G L 2 (]R d ), m G N Ubf^^U t C\ W™f(m d ,L 2 (Y d )), 

where, as usual, (x) = (1 + |x 2 |) 1//2 . 

The advantage of this construction is that the transformed Hamiltonian is a fibered operator 
over Y*. To assure that H = — A+Vr is self-adjoint in L 2 (IR d ) on the domain W 2 ' 2 (R d ), we 
make the following Kato-type assumption on the Pperiodic potential \RS\ Theorem XIII. 96]: 

(2.6) V v G tf oc (R d ) for d < 3, V r G L^ oc (R d ) with p > d/2 for d > 4. 
One checks that 

UbfHU b I = J dkH(k) 

with fiber operator 

(2.7) H(A;) = (_iV w + fc) 2 + *fr(i/), G R d , 

acting on the fc-independent domain T>q = W 2 ' 2 (T d ) C L 2 (T d ). Each fiber operator H{k) 
is self-adjoint, has compact resolvent and thus pure point spectrum accumulating at infinity. 
The eigenvalues are labeled increasingly, i.e. Eo(k) < E\(k) < E2(k) < . . ., and repeated 
according to their multiplicity. Since the fiber Hamiltonians are r-covariant, see |PST| . in the 
sense that 

(2.8) H(k + \) = t(X)- 1 H{k)T(\), VAeT*, 

the eigenvalues are T*-periodic, i.e. E n {k + A) = E n (k) for all A G P*. 

We denote by a*(k) the set {Ei(k) : n < i < n + m — 1}, k € Y* , corresponding to a 
physically relevant family of m Bloch bands, and we assume the following gap condition: 

(2.9) inf dist(<7*(fc), cr(H(k)) \ a*(k)) > 0. 

The following result borrows ideas in |C1H IC12"] and |Ne2| Theorem 3.3], where however a 
different definition of the Bloch-Floquet transform is used. 

Proposition 2.1. Let P*{k) 6 B(H{) be the spectral projector of H(k) corresponding to the 
set a*{k) C R. Assume that a* satisfies 12. 9\) . Then the family {P»(fc)} fceK d has the following 
properties: 

(Pi) the map k i— >■ P*(fc) is smooth from R d to B(Hf) (equipped with the operator norm); 
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(Pi) the map k i— > P*(k) extends to a B(Hf) -valued analytic function on the domain 

(2.10) n a = ^K£C d : |Im(Kj)|<a Vj G {1, . . . , d}| 

for some a > 0; 
(P2) the map k 1— > P*(k) is T-covariant, i.e. 

P^k + X) = t(X)~ 1 P^k)r(X) \/k G M. d , VA G T*; 

(P2) i/ie map k h- )• P*(k) is T-covariant, i.e. 

p*( K + x) = t(x)~ 1 p*{k)t(x) g n a , VAgT*; 

(P3) i/iere exists an antiunitary operato ^^\ C acting on T~L{ such that 
P*(-k) = C P*(k) C'- 1 and C 2 = 1. 

While properties (Pi) and (P2) are a consequence of the fact that the operator H commutes 
with the lattice translations, jointly with the gap condition (12. 9|) . property (P3) follows from 
the fact that the operator (II. ip is real, and corresponds to the time-reversal symmetry of 
the physical system. Notice that property (P3) generically does not hold true for periodic 
magnetic Schrodinger operators. 

Proof. One uses (I2.7P to define H(k) for every k G C d . Since for any kq G C d one has 

H (k) = H(k ) + 2(k - k ) ■ (-iVy) + (k 2 - k§)1 

and (— iV y ) is -ff (Ko)-bounded, the family {H(K,)} K£<C d is an entire analytic family of type (A) 
|RS[ Chapter XII]. Hence {H (n)} KeC d is an entire analytic family in the sense of Kato |RS[ 
Theorem XII. 9]. As a consequence, the set 

R = A) G C d x C : AG p(H(k))} 

is an open set, and (k, X) 1— > (H(k) — Al) _1 is an analytic function on R. Since R is open and 
(|2.9p holds, for every every ko G M. d there exist a positively oriented circle A(ko) C C and a 
neighborhood Afc C C d of feo, such that A(fco) separates <t*(k) from the rest of the spectrum 
and A(/co) C p{H(k)) for every k G Afc . Then the Riesz's formula 

(2.11) p 3t ( K ) = ^-(h (h(k) - z\y x dz 

defines a map k i— > P*(k) which is analytic in Afc . For k = k G Afc n it agrees with the 
spectral projection P*(k), so (Pi) holds true. 

Since cr* is r*-periodic, the circle in (|2.1ip can be chosen so that A(fco) = A(&o + A) for 
every A G T*. Thus, property (P 2 ) follows from dSHJ) and (f2~TTj) . 

Formula (|2.1ip yields a i3(%f)-valued analytic function on Ufc Afc D K d . In view of (P2) 
and the unique continuation principle, one may assume A/ co+ ^ = Afc for every A G T*, so 
by the compactness of Y* there exist a > such that U^A/^ D f2 Q , i.e. (Pi) holds true. 
Moreover, in view of (P2) and the unique continuation principle, one gets (P2)- 

W By antiunitary operator we mean a surjective antilinear operator C : H — > H, such that (C(p, Ci/j) H — 
(ip,<p) n for any (p,tp € H. 
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Property (P3) corresponds to time-reversal symmetry. This symmetry is realized in L 2 (M. d ) 
by the complex conjugation operator, acting as (Tip)(x) = ip{x) for tp E L 2 (M. d ). One directly 
checks that f = U BF TU^l acts on U T ~ L 2 (Y* ,%{) as 

(f<p)(k)=C<p(-k), L 2 (Y*,nf), 

where C is the complex conjugation operator in H{. From the fact that H commutes with T, 
taking into account ()2.9|) . property (P3) follows, see |Pa| for details. □ 

2.2. The Wannier functions and their localization properties. 

Case I. Simple Bloch band. We initially consider the case of a single isolated Bloch band, 
namely o~*(k) = {E n (k)} and fj2.9|) is satisfied, such that E n (k) is an eigenvalue of multiplicity 
one for every k. This case corresponds to the "simple direct isolated band " in |Ne2| . 

A (normalized) Bloch function corresponding to the n-th Bloch band is, by definition, 
any £ T-L T satisfying 

(2.12) #,-)GD H(k)ip(k, •) = E n (k)ip(k, •) and \\<p(k, -)\\ Hi = 1 VA: G Y*. 

Clearly, if ip is a Bloch function then ip, defined by tp(k,y) = e^^(p(k,y) for any function 
•& : — > R, is also a Bloch function. The latter invariance is often called Bloch gauge 
invariance. Notice that ()2.12p is equivalent to 

(2.13) P*(k)ip(k, ■) = <p{k, ■) and \\(p(k, -)||« f = 1 Vfc G Y*, 

where P*(k) is the rank-one projection on the eigenspace corresponding to the eigenvalue 
E n (k). 

Definition 1. The Wannier function w n G L 2 (M. d ) corresponding to a Bloch function 
<~p n £ % T for the Bloch band E n is the preimage of ip n with respect to the Bloch-Floquet 
transform, namely 

w n (x) := (UxF<Pn) (x) = — — r / dk e lk ' x ip n {k, [x]). 
v ' |y*|2 Jy* 

The translated Wannier functions are 

w n>y (x) := w n (x - 7) = T / dk e'^" 1 e lk ' x ip n {k,[x]), 7 G T. 

\Y*\2 Jy* 

Thus, in view of the orthogonality of the trigonometric polynomials and the fact that Ubf 
is an isometry, the functions {w ni7 } 7g r are mutually orthogonal in L 2 (M. d ). Moreover, the 
family {u; n , j7 } 7e r is a complete orthonormal basis of ti^p (Rani 3 *), where = fy* P*(k) dk 
is the total projector corresponding to the Bloch band E n . Notice that the previous definition 
and elementary properties do not rely on the gap condition (|2.9p . 

The (weak) localization of w n expressed by the fact that w n is in L 2 (R d ) is physically 
unsatisfactory, since it does not imply that the position operator and its powers have finite 
expectation value. Therefore, it is natural to search for conditions on (p n which guarantee a 
stronger decay of w n . In a nutshell, the localization properties of the Wannier function are 
determined by the regularity (continuity, smoothness, analyticity, . . . ) of the corresponding 
Bloch function. By a simple integration by parts, and an exact cancelation of the opposite 
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boundary terms, one obtains the following lemma. Here Xj, for j G {1, . . . , d}, denotes the 
the j th component of the position operator, i.e. (Xjip)(x) = Xj ip(x) for ^ in a suitable dense 
subspace of L 2 (R d ) 

Lemma 2.2. Let (p G TL, T , <p G C 1 (M. d ,'Hf) and let w =U~^ip. Then Xjw is m L 2 (R d ) and 

UsFi-iXjw) = d kj ip. 

By iterating the previous lemma, and taking into account (|2.5p . one concludes that if <p G % T 
is in C°°(lR d , %{), then the corresponding Wannier function decreases faster than the inverse of 
any polynomial, i.e. P(X\, . . . , Xd)w is in L 2 (M. d ) for any polynomial P. As for the exponential 
localization, by mimicking the proof of the usual Paley- Wiener theorem one gets the following 
result, see also |Kul| for a slightly different formulation. 

Proposition 2.3. Let VLp be defined by Ii2.10\) and 

n%/3 = {$ G I^JSlp,H{) : $>(z - A) = t(\)$(z) for all A G T*,z G . 
Lei 6e the restriction to M. d of a function $ G g analytic in the strip Qp. Assume that 

j \\<&{k + ih)\W dk < C Vh with \hj\ < f3 

with a constant C uniform in h. Then, the function w := U^p (ft satisfies 

[ e 2/3|x| \w(x)\ 2 dx < +oo. 

Sketch of the proof. Let k + ih be in Qg, and pose (j) h (k) := &(k + ih), so that 0^ G 7-L T . By 
shifting the integration contour by the method of residues, and by the r-equivariance of <p h , 
one gets 

w h (x) := (Wbf^) (x) = j dke ik - x $(Jfe + ih, [x]) = e h - x w(x). 
Then, by the unitarity of Ubf, 

<5>(k + ih)\\ 2 Hi dk < C. 

Since the latter constant does not depend on h for \hj\ < (3, one has He^'w; H^ngd) < +oo. □ 

Corollary 2.4. Let (j) be the restriction to M. d of a function $ G a analytic in the strip 
H a . Then, for every (3 < a, the function w := U^p <fi satisfies 

(2.14) f e 2/3 W \w(x)\ 2 dx < +oo. 

A function w G L 2 (U. d ) satisfying (|2.14p for some (5 > is said to be exponentially localized, 
while a function w G L 2 (R. d ) such that P(X\, . . . ,Xd)w is in L 2 (M. d ) for any polynomial P is 
said almost- exponentially localized. 

In view of the previous proposition, the existence of an exponentially (resp. almost-expo- 
nentially) localized Wannier function for the Bloch band E n is equivalent to the existence of 
an analytic (resp. smooth) Bloch function. Property (Pi) (resp. (Pi)) assures that there is 
a choice of the Bloch gauge such that the Bloch function is analytic (resp. smooth) around 




BLOCH BUNDLES, MARZARI-VANDERBILT FUNCTIONAL AND WANNIER FUNCTIONS 



9 



a given point. However, as several authors noticed |C11|, INe2] . there might be a topological 
obstruction to obtaining a global analytic (resp. smooth) Bloch function, in view of the com- 
petition between the regularity and the r-equivariance (remember that the Bloch function 
must be in H T by definition). This topological obstruction will be encoded in the concept of 
Bloch bundle, which we will introduce in the next subsection. Preliminarily, we describe the 
more realistic case of the composite Bloch bands. 

Case II. Composite Bloch bands. We consider now the generic case of a family cr* of m 
Bloch bands satisfying (|2.9p . Since the eigenvalues in the family <r* generically intersect each 
other, and the eigenprojectors of H(k) corresponding to single eigenvalues are not smooth at 
the intersection point, generically it is not even possible to find a system of locally smooth 
Bloch functions spanning RanP^(fc) at any k close to a given k$. Thus, the notion of Bloch 
function is relaxed and replaced by the following one |B1| IC11| . 

Definition 2. Let {P*(k)} kG ^d C £>(%f) be a family of orthogonal projectors satisfying (Pi) 
and (P2); with dimP^(A;) = m < +00. A function X £ 7~L T is called a quasi-Bloch function 
(for the family {P*(k)}) if 

(2.15) P*(k)x(k,-)= X (k,-) and X (k, ■) ^ V/c G Y*. 

A Bloch frame (for the family {P*(k)}) is a set {Xa}a=i,...,m of quasi-Bloch functions such 
that {xi(k), ■ ■ ■ ,Xm(k)} is an orthonormal basis o/RanP^A;) at (almost-) every k G Y* . 

A Bloch frame is fixed only up to a /c-dependent unitary matrix U(k) G U(m), i.e. if 
{Xa}«=i,...,m is a Bloch frame then the functions Xa(k) = YlbLi Xb(k)Ub, a (k) also define a 
Bloch frame. 

Definition 3. The composite Wannier functions corresponding to a Bloch frame {xa]^=i 
are the functions 

w a (x) := {U^pXaJ (x), ae{l,...,m}. 

As in the case of a single Bloch band, the existence of a system of exponentially localized 
(resp. almost-exponentially localized) composite Wannier function is equivalent to the exis- 
tence of an analytic (resp. smooth) Bloch frame. The topological obstruction to the existence 
of a regular (analytic, smooth or continuous) Bloch frame, already observed in |Clll INe2| is 
described in the next subsection. 

Remark 2.5 (Regularity of composite Wannier functions). We emphasize that, for Vr 
satisfying (I2.6p . the composite Wannier functions are actually in W 2,2 (M. d ) and, in general, 
one cannot expect better regularity properties. For example, if Vp has a Coulomb singularity, 
the Wannier functions are not smooth, as it happens for the eigenfunctions of the hydrogen 
atom. 

To study the regularity properties of w a , one notices that the corresponding quasi-Bloch 
function ip a (k) is in W 2 ' 2 (T d ) and 11^(^)11^2,2^) < C m ||V'a(^)||L 2 (T d )- Indeed, one can check 
that the previous inequality is a consequence of the fact that Vr corresponds to a multiplication 
operator in L 2 (T d ) which is A-bounded with relative bound zero, together with the fact that 
the union of the ranges of the functions {Ei : n<i<n + m— 1} is contained in a fixed 
compact set. As a consequence, in view of (|2.5p . one concludes that w a G W 2 ' 2 (W d ), namely 
w a is in the domain of H. 
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Remark 2.6 (Comparison with the physics literature). When comparing our definitions 
with the physics literature, one has to take into account that we are using a modified Bloch- 
Floquet transform, so that the fiber Hamiltonian (|2.7p has a /c-independent domain. This fact 
is convenient from the mathematical viewpoint. Alternatively, the classical Bloch-Floquet 
transform 

(2.16) (M(^):=^rE e " it7 ^+^ k G y G 

yields a decomposition Ubf {— A +Vr) = fj* H c \{k) dk where H c \(k) = — A -\-Vr acts in 
H k := E ^LO^) : + 7) = e ife "V(y) V 7 G T for a.e. y G R d } 

on the domain := Wj '<? (M. d )D'H] e . The unitary operator J = ZYbf^bf ac ^ s as (^^V) (k, y) = 
e lk ' y ip(k,y), is a fibered operator and its fiber, denoted by J(k), maps unitarily the space Hi 
into the space T~Lk- Since J{k)H c \{k)J{k)~ l = H(k), the operators H(k) and H c \{k) have the 
same spectrum {E n (k)} n ^. 

In the physics literature, usually a Bloch function is defined as an eigenfunction ip n (k, •) G 
T>k of H c \{k) = — A+Vr for the eigenvalue E n {k). By the unitary equivalence above, the 
so-called Bloch theorem, one has that tp n (k,y) = e lk ' y <p n (k,y), where tp n (k, •) is T-periodic and 
<p n G rl r . The latter is our Bloch function, as defined by (|2.12p . Consequently, in the physics 
literature the Wannier function is defined as w n = U^pijj n = Wgp(^ n , which coincides exactly 
with our definition. 



2.3. The Bloch bundle. To describe the geometric obstruction to the existence of a con- 
tinuous, smooth or analytic Bloch frame it is convenient to introduce the concept of Bloch 
bundle, following |Pa| . In this subsection we assume as given a family of orthogonal projectors 
{P*(k)}k£R d satisfying properties (Pi), (P2) and (P3), or their complex counterparts. Notice 
that we could also abstract from the specific case of the operator (|l.ip and consider such 
properties as convenient starting assumptions. Within this viewpoint, our approach can be 
applied to the periodic Pauli or Dirac operator, with obvious modifications. 

Proposition 2.7. To a family of orthogonal projectors {P*(fc)} fc6R d satisfying (Pi) and (P2) is 
canonically associated a Hermitian smooth vector bundle 8* over T* d , called the Bloch bundle. 
If (Pi) and (P2) are also satisfied, then is the restriction to = R rf /T* of a holomorphic 
Hermitian vector bundle £ it over Q a /T* . 

Proof. The idea of the construction is to firstly consider U fcgIR[ i Ran P*(fc) as a subbundle of 
the trivial bundle M. d x rl{ over ]R rf . Then we use the r-equivariance of the projectors to obtain 
a (quotient) vector bundle E* over the quotient space = R rf /r*. 

To construct one firstly introduces on the set M. d x Tif the equivalence relation ~ r , where 

(fc, if) ~ r (k', <p') <^ (k',<p') = (fc + A, r(A)~V) for some A G T*. 

The equivalence class with representative (k,(f) is denoted by [&,<£>]. Then the total space £* 
of the vector bundle is defined by 

£* := { [k, if] G {R d x Hi )/~ T : tp G Ran P, (k) } . 
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This definition does not depend on the representative in view of the covariance property 
(P2)- The projection to the base space ir : £ it — > is 7r[fe, <p] = u(k), where [i is the 
projection modulo T*. One checks that A is a smooth complex vector bundle with 
typical fiber C m . In particular, the local triviality follows from (Pi) and the use of the Kato- 
Nagy formula. Indeed, for any ko G M. d there exists a neighborhood Ok C M. d of kg such that 
— P*(ko)\\ < 1 f° r an Y k G Ofe . Then by setting (Kato-Nagy's formula |Kal Sec. 1.6.8]) 

(2.17) W(k) := (1 - (P*(k) - P*(k )) 2 y 1/2 (P*(k)P*(k ) + (1 - P,(fc))(l - P*(k ))) 

one gets a smooth map W : O h) -)• K(H { ) such that W{k) P*(k ) W(k)~ 1 = P*(k). If 
{Xa}a=l ...,m ls an y orthonormal basis spanning Rani^(/co) 5 then (p a (k) = W(k)xa is a smooth 
local orthonormal frame for £*, yielding the local triviality of the fibration A T*,. 

The vector bundle carries a natural Hermitian structure. Indeed, if V\,V2 G £* are 
elements of the fiber over fj>(k) G then, up to a choice of the representatives, one has 
vi = [k,(fi] and v 2 = [k, ^2], and one poses (ui,^)^ := (^1,^2)^. • 

As for the analytic case, one introduces an equivalence relation ~ T over £l a x Hf as above. 
Then the total space of £* is defined by 

£* = {[k, ip] G (ft a x % f )/~ r : ^ <= RanP„(K)} 

which is independent of the representative in view of (P2)- Local triviality follows again from 
formula (|2.17p . Indeed, for every Ko G Q a there exists a neighborhood N K0 C £l a such that 
||P*(/c) — -P*(fto)|| < 1 f° r every k G N Ko ; then, formula (12.17P yields a holomorphic map 
W : N K0 U(Ui) such that W (k) P*(k )W (k)^ 1 = P*(k). Notice that, since P*(k)* = P*(7c), 
the operator W(k) is not unitary when k ^M. d (1 N K0 . 

Since k t— > P*(k) extends k 1— >■ P*(/c) by (Pi), and both the maps are r-covariant, — >• 
is clearly a restriction of — > £l a /T* . The Hermitian structure over the vector bundle is 
defined as in the smooth case. □ 

The vector bundle £* is equipped with a natural u(m)-connection (Berry connection), in- 
duced by the trivial connection on the trivial vector bundle (M. d x Hf ) /~ r — > T^. The triviality 
of the latter is a consequence of the fact that U(T~Lf) is contractible whenever T~L{ is infinite- 
dimensional [Km]. 

The Bloch bundle encodes the geometrical obstruction to the existence of a global smooth 
(resp. analytic) Bloch frame. Indeed, the following result is implicit in |Pa] , 

Theorem 2.8. Let {P*(fc)} fcgK d C B(T~L{) be a family of orthogonal projectors satisfying prop- 
erties (Pi) (resp. (Pi)) and (P2) (resp. (P2)); with dimP*(/c) = m < +00. Then the following 
statements are equivalent: 

(A) existence of a regular Bloch frame: there exists a Bloch frame {Xa} a= \ m such 
that each Xa is i n C oc (M d ,Hf) (resp. each Xa is the restriction to M. d of a function 
Xa G T~(-T,a analytic on Q a ); 

(B) triviality of the Bloch bundle: the vector bundle £ if , associated to the family 
{P*(k)}] t£ ^d according to Proposition \2. 7\ is trivial in the category of smooth Her- 
mitian vector bundles over ( resp. the vector bundle £* is trivial in the category of 
holomorphic Hermitian vector bundles over f2 Q /T*). 
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The Bloch bundle can be trivial for reasons unrelated to time-reversal symmetry, as in some 
phases of the Haldane model |Ha[ IPr|. On the other hand, as a consequence of |Paj . under the 
assumption of time-reversal symmetry the Bloch bundle is always trivial in low dimension, as 
stated in the following result. 

Theorem 2.9. Let {P*(k)} k ^d C B(Hf) be a family of orthogonal projectors satisfying prop- 
erties (Pi), (P2) and (P3). Assume d < 3 and m > \, or d > 1 and m = 1. Then there exists 
a Bloch frame {Xa} a= \ m such that each Xa is real- analytic, i.e. Xa G C u (M. d ,7l{) flH T . 

Sketch of the proof . Since the family {P*(k)} k ^d satisfies properties (Pi), (P2) and (P3), 
in view of |Pa[ Theorem 1] there exists a smooth Bloch frame x = {xi)---iXm} C rl T H 
C°°{R d ,fli), i.e. the Bloch bundle is trivial as a smooth Hermitian vector bundle. Moreover, 
since (Pi) is also satisfied, {P*(k)} k ^d admits a holomorphic extension to £l a , see (|2.10p . 
Theorem 2 in |Pa| implies the existence of a family of holomorphic functions Xa '■ Hi, 
Xa G ?v a , such that {xi • • • ,Xm( K )} i s a (possibly non-orthonormal) basis of RanP*(/c) 
for every k G Q a . By restriction from Q a to M. d , and by a Gram-Schmidt orthonormalization 
procedure, one gets a real-analytic Bloch frame {xi, • • • , Xm} C C w (M d ,'Hi) fl rl T . □ 



3. The Marzari-Vanderbilt localization functional 



The long-lasting uncertainty about the existence of exponentially localized composite Wan- 
nier functions in three dimensions, settled only recently [BPCMj IPa| . and the need of an 
approach suitable for numerical simulations, lead the solid state physics community to ex- 
plore new paths. In an important paper |MaVa] . Marzari and Vanderbilt introduced the 
following concept. 

For a single-band normalized Wannier function w G L 2 (M, d ), one defines the localization 
functional by 

d r d / r N 2 

|2. 



w(x)\ dx 



(3.1) F MV (w) = VVar (Xj; \w{x)\ 2 dx) = \ \x^\w(x)\ 2 dx - V ( / xj 

which is well-defined at least whenever J" Kd |x| 2 |u;(x)| 2 d2; < +00. 

More generally, for a system of L 2 -normalized composite Wannier functions w = {w±, . . . , w m } C 
L 2 (M. d ) the Marzari-Vanderbilt localization functional is 

m m „ m d , „ \ 2 

(3.2) F M y{w) = y^-FMy(wq) = ^ / N 2 |w a (x)| 2 <ix - V] V] ( / Xj|w a (x)| 2 <ix ) . 

We emphasize that the above definition of Fmv(w) includes the crucial constraint that the 
corresponding Bloch functions (p a (k, •) = (Ubf w a )(k, •), for a G {1, . . . , m}, are a Bloch frame, 
i.e. {(pi(k, •),..., (f m {k, •)} is an orthonormal set in T~Lf for each G Y* and 

(3.3) Span c {tpi(k, ■),..., <p m (k, •)} = P*(k)(Hf) , Vfc G F* . 
According to Remark 12.51 the latter condition actually implies w a G W 2,2 (M. d ) = T>(H). 

Definition 4. Let {-P*(£0}fceR d ^~ &(flf) ^ e a family of projectors satisfying properties (Pi) 
and (P2), with dimP !f (/c) = m < +00. ^4 system of maximally localized composite Wan- 
nier functions is a global minimizer {wi, . . . ,w m } of the Marzari-Vanderbilt localization 
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functional -Fmv ^ n the set VV m := (T>(H) n T>(X)) m , under the constraint that {<fi, ■ . . , (p m }, 
for (p a = lA-Q^Wa, is a Bloch frame. 

A natural problem, raised in |MaVa| . is the following. 

Problem 3.1. Let {-P*(fc)}fceR d *~ ^(^f) De a family of projectors satisfying properties (Pi) 
and (P2)) with dimP„,(/c) = m < +00. 

(MVi) (Existence) prove that there exists a system of maximally localized composite Wan- 
nier functions; 

(MV2) (Localization) prove that any maximally localized composite Wannier function is 
exponentially localized, in the sense that there exists ft > such that fj2. 14j) holds. 

Since the (modified) Bloch-Floquet transform U^F '■ L 2 (M. d ) — > L 2 (Y*;rl{) is an isometry 
and it satisfies (UBFXjg)(k,y) = i^ur(UBFg)(k,y), the functional (13.21) can be rewritten in 
terms of the Bloch frame (p = {<p±, ■ ■ ■ , as 
(3.4) 

Correspondingly, in view of (J23J), the space W = V(H) n V(X) = V(H) n T>({X)) is mapped 
by the Bloch-Floquet transform into 

n T n Lf oc (R d , w 2 ' 2 (T d )) n w^(R d , L 2 (r d )) =■. w. 



Hereafter, to solve problem (MVi) for any d, we will make the following 

Assumption 1: there exists a Bloch frame \ = {Xli • • • ; Xm} C rl T such that Xa S W. 

In the case most relevant to us, i.e. Schrodinger operators for d < 3, the previous assumption is 
automatically satisfied, since Proposition 12.11 and Theorem 1 2 . 9 1 provide the existence of a Bloch 
frame which is even real-analytic. While this extra regularity is unessential for problem (MVi), 
it will be crucial when dealing with (MV2). Notice that, under the previous assumption, the 
set of admissible Bloch frames in Definition 0] is non-void, so problem (MVi) makes sense. 

By using the previous Bloch frame, we can lift the functional (|3.4p to W 1,2 -maps from 
to the unitary group U(rn), i.e. to T*-periodic maps from M. d to 14 (m). Indeed, given 
any map U E W ' (Tj,W(m)) one defines a Bloch frame ip = {tp\, . . . , tpm} C VV by setting 
tp = X-U, i.e. ip a (k, •) = Yjb Xb(k, •) U b>a (k). Vice versa, if if = {ip 1 , ... , ip m } C VV is a Bloch 
frame, then pointwise <p a (k,-) = Y.bXb{k,-)U b , a (k) with U b>a {k) = (xb{k) , f a (k)}, hence 



dy 



dk 



^Pa{k,y)i^ L 
T d okj 



(k, y) dy 
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For the given reference frame x, the functional (|3.4p in terms of the gauge U becomes 



(3.5) 



(U; X ) 



d 

E 

d 

E 



tr 



tr 



dU* 



+ m 

0=1 



dXa(k,-) 


2 


dkj 





life + 



5C7 
9fei 




Here the matrix coefficients Aj G L 2 (T^;u(m)) are given by the formula 



dk 



(3.6) 



[A,-(AO] 



cb 



Xc(k, 



dXbjk, ■ 
dkj 



When x is real- analytic, the functions Aj G C w (T^;u(m)) represent the antihermitian connec- 
tion 1-form induced on the (sub)bundle by the trivial connection on the bundle M. d xT~Lf- 

Moreover, 
(3.7) 

{wi, . . .,w m } C W 



inf < FmvH 



Ubf w is a Bloch frame 



mi{F MV (U; X ) : U G W 1 ' 2 ^; U{m))} . 



Therefore, problem (MVi) is equivalent to showing that the r.h.s. of (|3.7p is attained. Anal- 
ogously, in view of Corollary 12.41 problem (MV2 ) corresponds to show that any minimizer of 
Fmv{''-,x) i s real- analytic, provided that x ls a l so real-analytic. 

Remark 3.2 (Rough reference frames, as in numerical simulations). When reformu- 
lating problem (MVi) in terms of the functional Fmv{U',x) the regularity of the reference 
frame x plays no essential role, provided x is in W. Indeed, in view of (|3.7p . the infimum on the 
r.h.s does not depend on the choice of X- Moreover, if for a particular choice of x the infimum 
is attained at U, then for another choice x C W the infimum is attained at U = V U, where 
X = X ■ V. The matrix V, defined by Vb >a (k) = (Xb(k) , Xa(k)), is in W 1 ' 2 (T^;U (m)) , hence 
U is also in W 1,2 (T^;U(m)). This observation justifies the fact that a minimizer (whenever it 
exists) can be evaluated starting from any reference Bloch frame x C W, even a discontinuous 
one, as it happens in numerical simulations. 

Remark 3.3 (Minimizing over smooth gauges). To compute the infimum of Fmv{U',x) 
it is sufficient to consider smooth change of gauges. More precisely, for d < 3 and for any 
fixed Bloch frame x C W one has 

inf [P MV (U;x) : U G W 1 ' 2 ^; U{m))) = inf {f mv (U; X ) ■ U G C°°(TS;^(m))} . 

As a consequence, in numerical implementations, to compute the above infimum one can let 
U vary in any set S such that C°°(T*;^(m)) C S C W 1 ' 2 ^; W(m)). 

This result follows from the strong density of smooth maps C ao {T* d ;lA{m)) C W 1,2 (T^;U(m)) 
for d < 3. If d = 2 the latter claim is essentially the approximation by convolution followed by 
the nearest-point projection onto the unitary group; if d = 3, the claim follows from the fact 
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that the homotopy group ^(^(m)) is trivial and from the fact that T3 has the 1-extension 
property with respect to U(m), see |HL[ Theorem 1.3 and Section 5]. 



3.1. Existence of minimizers. The following results shows that the right hand side in (|3.7p 
is attained. The proof is a simple modification of the direct method in the calculus of varia- 
tions, in order to handle a natural invariance of the functional (|3.2p . Indeed, if {w±, . . . ,w m } 
are composite Wannier functions satisfying f)3.3[) and {71, . . . ,j m } C T, then 

(3.8) F M v(wi, ■ ■ ■ ,w m ) = F M v(wt, . . . ,w m ) , w a (x) = w a (x + -f a ) , 1 < a < m . 

Moving to Bloch functions, we have (p a {k, •) = (Ubf Wa^j (k, •) = e lk " <a (Ubf 'Wa^j (k, ■), so that 

{ipi(k, •),..., <f m (k, •)} is still orthonormal in T~Lf and ()3.3|) holds. Correspondingly, the func- 
tional (|3.5|) has the invariance 

(3.9) F MV (U; x) = F MV 0; x) , U(k) = diag (>>) U(k) . 

Theorem 3.4. Let {P*(k)} kG ^d C £>(%f) be a family of orthogonal projectors satisfying prop- 
erties (Pi) and (P2); with dimP*(A;) = m. Assume that there exists a Bloch frame x C VV. 
Then there exists U E W 1,2 (T* d \U{m)) which is a minimizer on W ' (T^;W(m)) of the local- 
ization functional Fmv (' j x) defined by \3. 5\) . 

Proof. Note that F(-;x) is nonnegative, because of (13 .7h and (13. 2p . Note also that E(U) = 

X^=i fj* \ tr ( 7^(^)mT"(^) ) i s the standard Dirichlet integral so it is (sequentially) weakly 

lower semicontinuous on W 1,2 . Thus F(-;x) 1S weakly lower semicontinuous in W 1,2 , because 
the first term in (|3,5p is, up to a constant factor, E(-) and the other terms are clearly weakly 
continuous (because of the compact embedding W 1,2 (T^) =— )■ L 2 (T^)). In order to apply 
the direct method it remains to show that the functional is coercive, so that there exists a 
minimizing sequence which is bounded in W 1,2 (T^;U(m)). However, coercivity clearly fails 
because of the natural invariance of F(-; x) given by (|3.9p (the W^^-norm of U can be made 
arbitrarily large as \^ a \ — > 00). 

In order to fix the argument it is enough to take advantage of this invariance in the form 
(|3.8p . Indeed, for any admissible U we can choose 71, . . . ,7m so that the corresponding U, 
defined as in (I3.9p . gives the same value of Fmv ("j x) but in such a way that the corresponding 
Wannier functions {w a } satisfy Y2j=i \Sr<i ^jl^a^)! 2 ^) < C for some absolute constant C > 
independent of U. Thus, if along a sequence {U n } the functional Fmv is bounded, up to a 
suitable choice of the translation parameters {7"} we have a uniform bound for the modified 
sequence {U n } in W 1,2 (the third line in (13.51) is now bounded by C 2 md and the first two are 
easily seen to be equivalent to the Dirichlet energy up to additive and multiplicative constants 
depending only on A, C and m). 

Then, up to subsequences, U n weakly converges to some U S W 1,2 {T* d ;lA{m)) which is a 
minimizer. □ 
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3.2. The Euler-Lagrange equations. In this subsection we derive the Euler-Lagrange equa- 
tions corresponding to the functional (I3.5p . To our knowledge, these equations are new in the 
literature. We consider the unitary group U(m) as isometrically embedded into the space 
M m (C) with the standard real Euclidean product (A , B) = Retr(A* B). Recall that U{m) 
is a compact real Lie group, the induced metric is biinvariant, its Lie algebra is given by the 
real vector space of complex antihermitian matrices and it has dimension m 2 . 

Let ip e C°°(T* d ; M m (C)) and for e / fixed let U{k) + £<p(k) be a free variation of U in the 
direction if. In a sufficiently small tubular neiborhood O of U{m) in M m (C) there is a well 
defined smooth nearest point projection map II : O — > U(m), so we can consider the induced 
variations 



(3.10) U £ {k) := U(U(k) + eip(k)) = U{k) ( I + e- [U^(k)ip(k) - {U~ x {k)^p{k))*} ) + o{e) 



Simple calculations on each term in (13.5() yield 



(3.11) 
d 
de 



tr 



=o^t 



(3.12) 

and 
(3.13) 



1 ' dkj 

d_ 

de 



tr 



/ tr 


¥>*(*){ 


It* 





dU | 
dkj 



Jr* 


tr 


_dk^ 


dU* 

dkj 


(k) 


dk 3 



dU 

dkj' 



dU 
dkj 



dk 



dU 
dkj 



dk . 



dA ■ dU 



dk . 



d 

ds 




dU 

[^±(k) + A 3 (k)U £ (k) 



dk 



tr 



|^(fc) + Aj(k)U(k)J <P + U(k)GPU-\k) (^-(k) + Aj(k)U(k)j } 



Here the constant (purely imaginary) diagonal matrices {G 3 } C M m (C) are defined as G 3 
diag (/ T . U*(k) \§jjL{k) + Aj(k)U(k)] dk), where [diagM] o6 = M aa 5 ab . 



dk 



Thus a map U 6 W 1,2 (Y* d ]lA{m)) satisfies —Fmv(U e ; x)|e=0 = if and only if U is a weak 



solution to the Euler-Lagrange equations 



(3.14) 



j-i j j j =1 

d 



d^ A]U dk 3 u A 'dk, 



+ 



( dU 



- 



\dkj 



+ AjU G J + UG J U~ 



i f du 

dkj 



+ AjU 



0. 
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4. Continuity of minimizers 



The goal of this section is to show that a change of gauge U G W 1,2 (T d ;U(m)) which 
minimizes the functional (|3.5p is continuous, whenever \ ls real-analytic. This result will be 
true if m = 1 for any d, if 1 < d < 2 for any m and for d = 3 under the restriction 2 < m < 3. 
The proofs in the three cases are different. The first case, treated in the next proposition, is 
the simplest and it gives even real-analiticity. In order to deal with the other two cases, we 
will need a dimension dependent argument which will occupy the rest of the section. Note 
that in the first two cases continuity holds for any solution of the Euler Lagrange equations. 
In contrast, in the three dimensional case continuity relies on energy minimality in an essential 
way. 

To deal with the regularity of the minimizers, we make the following assumption which, as 
already noticed after stating the weaker Assumption 1, is automatically satisfied for d < 3 or 
m = 1 (compare Theorem 12.9ft . 

Assumption 2: there exists a Bloch frame x = {xi>---)Xm} C 7-L T such that \ a G 
C u (R d ,H t ) for every a G {1, ... ,m}. 

Proposition 4.1. Let d = 1 and m>l, or d > 1 and m = 1. Let U G W 1,2 (T* d \U{m)) be a 
weak solution to equation (|3.14ft . Then U is real- analytic. 



Proof. Assume d = 1. Since U G W 1 ' 2 ^, U(m)), then (l3~14jl yields j^U G L 1 . Thus U G 
W 2,1 . Therefore U is in C 1 by Sobolev embedding, and (13.14ft implies U G C 2 . Analogously, 
if U G C n+1 solves (I3T41) . n > 1, then U G C n+2 , so by induction U e C°°. Finally, since 
A,- G C w for all j, by the standard ODE regularity theory one obtains U G C w . 

Now assume m= 1. Since W(l) is abelian the equation (|3.14ft reduces to 




A dU i at/ / A OA 

(41 » -^+E^f-'^= Ey* 

j=i j j \j=i 

Recall that in any sufficiently small ball BcTJwe have U(k) = e 1 ?^ for some / G W 1,2 (B; R) 
(see |BZ| ). Thus, equation (|4.1ft reads A/ = i X)^=i §^ G C""(JB;R), because the matrices 
in the Berry connection are antihermitian and real- analytic. Since the Laplacian is analytic- 
hypoelliptic we conclude / G C W {B) and in turn U G C W (T^;ZY(1)) since the ball B can be 
choosen arbitrarily. □ 



4.1. Continuity in the two dimensional case. We are going to prove continuity of any 
weak solutions to (|3.14ft in the case d = 2. The argument here is just sketched, since, up 
to a standard localization argument, it essentially follows the proof of continuity for weakly 
harmonic maps from a two dimensional domain into spheres (see e.g. |LW1| . Chapter III, 
Section 3.2 pag. 57-61). 

We start with the following auxiliary result which is a straightforward consequence of 
jCLMSj . 
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Lemma 4.2. Let d > 2, m > 2 and let U G W l > 2 {T* d ; U{m)). Assume W G L 2 (T*;u(m)) for 
j G {1, . . . ,d} and div B = in T>'(T* d ). If U is a weak solution to 

d 



(4.2) AU = + f and U~ x f G LP(T*; u(m)) 



/or some p > 1, i/ien f7 G VF 2,1 (T^; U{m)). 
Proof. We have 

m <i ^g^y- m 

6=1 j=l 3 b=l 

where for fixed a,b,c the vector fields B^ and = VU a b are in I? and satisfy div7?(, c = 
and curl.E a 6 = in the sense of distributions. According to [CLMSj, we have E a ^ ■ B^ c G 
^iocC^d) and Uc G "Wioc(Trf), i.e. the right hand side of (02]) is in the local Hardy space 
^loc( T S) c Ll ( T *d)- Here 9 G ^loc( T d) means tliat 9 G U\ oc {B) for every sufficiently small 
ball i? C (namely, for every ball with radius smaller than the injectivity radius of the 
exponential map). For the definition and the basic properties of HL(i <i ) we refer to [St| 
Chapter 3]. Thus the conclusion follows from |LWll Theorem 3.2.4]. □ 

Based on the previous lemma we have the following intermediate result. 
Proposition 4.3. Let d > 2, m > 2 and U G W 1 ' 2 ^; U(m)). If U is a weak solution to 

(4.3) AU = ^ U ' 1 ^. + / and U ~ X f G L 2 (T* d ;u(m)), 

j=i i i 

then U G W 2 ^(T* d ;U(m)). In particular, if d = 2 then U G C°(T* 2 ;U(m)). 

1 ( tt* dU dU* ; 



Proof. If we set B 3 = ± (U*§£; - %jU\ then B> G L 2 (T*; u(m)) for j G {1, . . . ,d}, U is 

a weak solution to gj2]) but div 5 = §(17*/ - fU*) ^ in V(T* d ). Note that div B G 
Lq(T^; u(m)), the space of zero-mean L 2 -integrable functions, hence A" 1 div 7? G W 2,2 (T d ; u(m)) 
by elliptic regularity and if we set 

B = B- VA -1 divS, / = / + VC/ • VA" 1 divS , 

then 7?, / and C7 satisfy the assumptions of Lemma [4.21 by Sobolev embedding (for some p > 1 
depending only on d), hence we obtain U G W 2,1 (T d ;U(m)). Finally, if d = 2 the improved 
Sobolev embeddings into Lorentz spaces yield V£7 G 7 2 ' 1 and in turn U £ C° (see |LW1| . 
Theorem 3.2.7 and 3.2.8 respectively). □ 

Going back to equation (|3,14p we have the following important consequence 

Corollary 4.4. Let d = 2, m > 2 and U G ,2 (T*>;W(m)) a weak solution to equation 
(j33H) . Then U G W 2 ' 1 ^; U{m)) and U G C°(T|;W(m)). 
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Proof. If U G W l > 2 (Y* 2 ;U(m)) is a weak solution to equation (|3.14|) then U satisfies also 
for a suitable / depending on U . Here t/" 1 / G L 2 (Y 2 ] u(m)) because of the regularity property 
of £/ and the Berry connection A,- G C^ClT^, u(m)), in view of the structure of equation (|3.14p . 
Thus, the conclusion follows from Proposition 14.31 □ 

4.2. Continuity in the three dimensional case. The goal of this subsection is to show that 
if d = 3 then minimizers of the localization functional (j3.5|) are continuous, at least if m < 3. 
Roughly speaking the idea is to prove that at smaller and smaller scales minimizers look like 
minimizing harmonic maps into the unitary group U{m) which are degree-zero homogeneous. 
Since the latter are constant, at least for m < 3 (see Corollary I A. 101 below) . then the former 
are continuous at a sufficiently small scale. All the techniques in this section are inspired by 
the regularity theory for minimizing harmonic maps (see |SU1| . [ST] and |LW1| ). 

The first condition we need to study minimizers at small scales is the stationarity condition 
with respect to inner variations. 

Lemma 4.5. Let d > 3 and let U G W l,2 (Y* d ; U{m)) be a minimizer of (13. 5ft . For each 1 < 
j < d define constant diagonal matrices G j = diag f t/*(J^- + AjU)dk. Let $ G C°°(T* d ;R d ) 
be a smooth vector field and *$> e (k) = k + £&(k) be a family of diffeomorphisms (for e small 
enough). Then 



d 
de 



e=0 



F MV (Uom £ ; X ) =J2 



:)■<■ 



d<$> c 

dk r . 



tr 



dU* dU 



dU 



S k 1 3k j +2 9k 1 UA '- 2(PU ^ + A > U 



dk- tT 



.)■<■ 



dU* dU dU* dU 
+ 



(4.4) 



dk c dkj dkj dk c dk 



2^U*A j + 2G j U*^- 

Okr. 



(dU_ 
[dkj 

+ 



+ 



.M' 



' dU dAj 
dkj dk r 



dk a 



Proof. First note that def-D^ 1 (k) = 1 — ediv$ + o(e) uniformly on T%. As a consequence 
the change of variable and the chain rule give 



(4.5) 



d 
ds 



e=0 



f(^ e (k))g(k)dk 

for any / G L 1 and for any g G C . Now we set U £ (k) = U o ty £ (k), so that 
(4.6) 



f(k)g(k)dw^(k)dk-J2 I f{k)^{kW(k)dk 



dk c 



dk 



Thus, a repeated application of (14. 5|) and (14. 6|) gives 



(4.7) 



£ = 



E 



3,c 



dU* dU £ 
dkj dkj 



d$ c (k) 
dkj 



dk 



^ j div®(k) tr 



dkj 



dU* dU £ 
dkj dkj 



+ 



tr 



Of/* dU dU* dU 
dk c dkj dkj dk c 
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(4.8) 



d 

ds 



e=0 



£/ te 



3 

dU* dU 



dkj dkj " 



dk 



U 



dkj dkj 



. n .d$ c dAj, c 



+£/* 



<9/c r . dk c 
d<5> c 



dk+ 



dkj 



dk. 



Similarly, taking the definition of G 3 into account, one has 

'du, 

£=0 



(4.9) 



d 

de 



£ 

j,a 



v 3 



3,c 



dk c ^ ^ dk c dkj G'U q^J^ 



Combining (14. 7p . (I4.8P and (14. 9|) with the definition of (13. 5p and reorganizing the sum we 
easily have (j4~4|) . □ 



The first consequence we obtain is a sort of perturbed monotonicity formula (in the spirit 
of the monotonicity formula for almost harmonic maps; see |Mos] . Chapter 4). 

Proposition 4.6. Let d > 3 and let U G W 1,2 {T* d ;U{m)) be a minimizer of the localization 
functional (13. 5ft . T/ien there exist C > and R > such that for each ko G and /or eac/i 

i?2 wii/i < Ri < R2 < R we have 
(4.10) 



1^2 / l VC/ | 2 + / 

JB Rl (k ) JR^k-ko 



<R 2 



\k - k \ d - 2 



dU 



dr 



< 2CR 2 + {l + 2CR 2 )-^ [ \VU\ 2 . 

R% JB R2 {k ) 



As a consequence, the function R 2 d 
is finite and vanishes as R — > + . 



|VJ7| 2 dfe has a limit and ^B R (ko) \^rU\ 2 \k— ko\ 2 d dk 



Proof. Clearly by translation invariance it suffices to consider the case ko = 0. We take 
h G C°°(R) an increasing function such that hit) = for t < and h(t) = 1 for t > 1. If we 
take $(/c) = kh(R- \k\) in (JO) we have div$(A;) = dh(R- \k\) - \k\h'(R- \k\) (d G N is 
the dimension) and 



(4.11) 



d ® c _ rei. fc c^ , dU d$ c _dU dU , ^ 9ij 
d% J "jjfep ' ^dk c ~dk~ ~ dk~ k i~dr~ h ^^ 



dk r 



dA, 
1 dr 



h. 



Combining (|4.4p and (14. lip we get 
= ^2 [ {dh- \k\ti)tr 



dkj dkj dkj 3 



dU 
dk,-, 



+ A S U 



£ 



fcj tr 



dU^dU_ + dU^dU_ 2 9U_ u * A + 2G i u *?E 
dr dkj dkj dr dr 3 dr 



h'+ 
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(4.12) + 2^ / tr 



dU*8U dU TTifA ^ jTT *dU ...dUdAj n ., TT *dA jT ' 



dkj dkj dkj " 



dkj dr 



Or 



h . 



Passing to the limit as h — > X{t>o}i as m (HE Chapter 2], yields 



(4.13) = [ tr 

i J B R 



(2-d) 



OU* dU 
dkj dkj 



tr 



dB f 



dU* dU 
dkj dkj 



2R I 

ldB R 



dU* dU 
dr dr 



dB R 



\k\ti 



+ 



2 E/ 1 

'Br 



(d _ i) (^U*Aj - GW*^f] - dXPXTAjU + |fc||^^_T " \k\G^U*^U 



\dkj 



dki 



dkj dr 



Oi 



Note that R d ~ 

for a.e. R > 0. Thus, dividing in (|4.13p by and integrating on [R\,R 2 ] we have 



dR \ R2 d Ib r / J = ( 2 ~ d ) lB R f- R IdB R f for an y hitegrable function / and 

?<f-l 



(4.14) 



1 



nd-2 . 

K l Jb Ri 



\vur + 



Ri<\k\<R 2 



Id— 2 



dU 



Or 



1 



i?: 



cf-2 



|vc/p 



R 2 



JR 1 <\k\<R 2 K J B] 



W(k,U), 



where 

V(k,U) = 2\ 
and 

W{k,U) = 2tr 



tr 



dU TT>k . ^.nrr* ( 9U . TT \ kj dU . 



_ _2 Gijj* — 
dr 



\k\ 



dU dU 
(d-l)i ^-U*Ai - GPU* 



dk^ 



dk. 



ftTT ft A ■ ft A ■ 

dGW*A j U +\k\^-^P- - \k\GW*^U 
dkj dr dr 



Notice that the first line in (|4.14p is the usual monotonicity identity for harmonic maps from 
with values into U(m). To estimate the extra terms in (I4.14p . observe that for fixed matrices 
Cr 7 we have \V(k, U)\ < C\k\ 2 ~ d (l + | V£/|) where the constant depends on the C° norm of the 
Berry connection matrices Aj. Similarly we have |W(fc, U)\ < C(l + |V£/|) where the constant 
depends only on the C 1 norm of the Berry connection matrices. Thus, if Iq > 1 is an integer 
and Ri > 2~ l °R 2 , R% < 1 and < 5 = R 2 < 1 we have (writing f B for r^r f B ) 



[ \V(k,U)\<C [ |A;| 2 - d (l + |VC/|)< 

JR 1 <\k\<R 2 JR 1 <\k\<R 2 



iO-1 n 

CRI + CY, / 

1=0 " 1 



1 R 2 <|fc|<2- ! i? 2 



io-i r~n 

vu\<CR 2 2 + cY, 2 ~ 2l RlJf |VC/| 2 < 

1=0 V 2-lRn 



CR\ + 5 sup —j-^ 

re[2- l 0R 2 ,R 2 ] r J 



\VU\ 2 + C5~ l Rl < CR 2 1 + 



1 

€[2- l 0R 2 ,R 2 ] r JB r 



\vu\ 
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On the other hand, a similar estimate with 5 = Si = 2~ l R 2 gives 
" R2 dR f „„„ * § f R2 dR 



\W(k,U)\<CR 2 + \ 

JR 



i?i K Jb r Jr.! -K 7s h 



VJ7I < Ci?^+ 



'0-1 r 2-'i? 2 jo i /■ 



Jo-l 



'Zo-1 



\Z=0 / rG[2- ! 0R 2 ,J?2] ' ■'■Br 



212-^2 
l2-'- 1 -R 2 




< Ci? 2 1 + SUp 

Combining the two estimates above we obtain 
(4.15) 

f \V(k,U)\+ [ R2 -^ [ \W(k,U)\<CR 2 (l+ sup *f |V/7| 2 ] 

JRK\k\<R 2 JRi n JB R \ re[2- l 0R2,R 2 ] r •> B r J 

Going back to (|4.14p . if R 2 < R ■ = min{l, ^C^ 1 } is fixed, taking the supremum over R\ £ 
[2~ l °R 2 , R2] and estimating the right hand side using (|4.15[) one easily obtains 

fP i/ |VC7| 2 < 1 + 2— / |VC/| 2 . 

Combining (|4.14p . (|4.15|) and the previous inequality we finally obtain (|4.1U|) . Letting R\ — > 
in (j4TTU|) we see that R 2 ~ d f B \VU\ 2 dk is bounded and f B \d r U\ 2 \ k\ 2 ~ d dk is finite and 
therefore vanishing as R — > + . As a consequence, letting R\ — > and R 2 — > in (|4,10p . it is 
straightforward to see that R 2 ~ d J Br \X7U\ 2 dk has a limit as R -> 0. □ 

Remark 4.7. A simple consequence of (|4.10p is that if lim^o S 2_rf fB R (k ) \W\ 2 dk = e 

then there exists Rq > such that su Pfc6Bfl (feo) su Po<-R<_R -^- 2d lB R (k) \^U\ 2 dk < 2e, i.e. 
at sufficiently small scales the scaled energy is locally uniformly bounded by its limit at any 
point. 

The second ingredient is a compactness theorem for the scaled maps Un{k) = U(ko + Rk) 
which is similar to the compactness theorem for minimizing harmonic maps (see |SU1| and 
[51], Chapter 2). 

Proposition 4.8. Let d > 3 and U £ W l ' 2 (T*,; U(mj) be a minimizer of the localization 
functional (|3.5|) . If we define Un{k) = U(ko + Rk), R > 0, then up to subsequences Ur — > Uq 
strongly in W lo ' c (M. d ;U(m)). In addition Uq is a locally minimizing harmonic map and it is 
degree-zero homogeneous. 

To prove the previous result we need the following two auxiliary lemmas. The first is a 
simple consequence of local minimality. The second is a nonlinear interpolation lemma due to 
Luckhaus which we state in our specific context (for the general case see e.g. |LW1| . Lemma 
2.2.9; see also [Sj], Chapter 2 for a proof). 

Lemma 4.9. Suppose U and Ur as in Proposition \4-8\ Let p G (0,1) and let {ur} C 
W 1,2 (B p ;U(m)) a bounded sequence such that Ur = vr on dB p for each R > 0. Then 



liminf/ \VU R \ 2 < liminf / IV 
R^o Jb J Bo 



vr\ 2 
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Proof. Define vn(k) = vji(R~ 1 (k — k )) so that vr G W 1,2 (B p ji(ko);U (m)) . Since vr = U on 
dBpft(ko) we can extend them as U to the whole T^. The assumption on vr and Proposition 

1J] clearly imply / \VU\ 2 + \Vv R \ 2 = 0{R d ~ 2 ) as R ->■ 0, hence formula (J53|) and 

JB pR (k ) 

simple calculations using Cauchy-Schwartz inequality give 

(4.16) F mv (v r ;x)-Fmv(U; X )= [ \^v R \ 2 - f \X7U\ 2 + 0(R d - r ) , 

■>B pR (k ) JB pR {k ) 

because vr and U coincide outside B p R(ko). Since U is a minimizer of Fmv the right hand 
side of (|4.16|) is nonnegative, hence scaling back and taking the definition of Ur and vr into 
account the conclusion follows as R — > 0. □ 

Lemma 4.10 (Luckhaus). Let d > 3, m > 2 and let u, v G W 1 ' 2 (5 <i-1 ; U(m)). Then, for 
each A G (0,1) i/iere is w £ W 1 ' 2 (S d ~ 1 x (1 - A, 1); M m (C)) such that w\ S d-i x ^ = u, 

Hs<*-i x {i-A} = v > 

(4.17) f \Vw\ 2 <Cx[ (\V T u\ 2 + \V T v\ 2 ) +CX- 1 [ \u-v\ 2 

,/,S<*-ix(l-A,l) Js^ 1 J S 1 *- 1 

and 

i i 

dist 2 (w()t),W(?n)) < CA 1 ^ H (|V T u| 2 + \V T v\ 2 )^j * (f \u - v\ 2 ^j % + 

(4.18) +C\- d [ \u-v\ 2 

for a.e. k G S^ 1 x (1 — A, 1). i/ere is the gradient on S"^ -1 . 



Proof of Proposition ^ -8\ We essentially follow the proof of |LW1| . Lemma 2.2.13, with mi- 
nor modifications. By Proposition 14.61 up to subsequences we may assume Ur — Uq in 
W 1,2 (B\;U{m)) where Uq is a degree-zero homogeneous map. Thus, it is enough to show 
strong convergence and minimality in some ball B p C B\ to get the same properties on any 
Bp C W 1 for any p > 0, by scale invariance of Uq and the existence of the full limit of 
R 2 ~ d J Br \VU\ 2 asi?^0. 

Let B\ C M. d and 5 G (0, 1) a fixed number and let w G W 1 ' 2 (Bi;U(m)) such that w = Uq 
a.e. on £?i \ Bis- By Fatou's lemma and Fubini's theorem, there is p G (1 — 5, 1) such that 

lim f {UR-UofdH'*- 1 = 0, / (|VC/ R J 2 + |W | 2 ) d-H^" 1 < C < oo . 

JdB p JdB p 

Applying Lemma [4.101 to A = Xr < 5, u = Ur(p-) and v = w(p-) = Uo(p-) for a de- 
creasing sequence of numbers Xr — > 0, we conclude that there exists a sequence of maps 

/ \ i/2d 

w R G W 1 ' 2 (B p ;M m (C)) such that if we chose e.g. X R = I f \U R - UofdH^ 1 ) < 5, 



*dB p 



'p 

then we have 



w R (k) 



U R (k), \k\=p, 
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Afl / {\V T U R \ 2 + \V T U \ 2 ) + A^ 1 f \U 

JdBp JdBp 



(4.19) 

/ iVwjil 2 < C 

and dist(u>_R,W(m)) — >■ uniformly on B p \ B p ^i_x R ) as i? 
{v R } C W^B^M) by 

(4.20) v R (k) 



U n \ 2 



i?->0 



0. Define comparison maps 



w 



<p(l-Afl), 



nO fl (fc)), - Ajz) < |A;| < p, 

where II : O — > U{m) is the nearest point projection. Then, by minimality of Ur, Lemma [479 
and fj47T9|) - (f4~20|) we obtain 

[ |W | 2 < liminf f \VU R \ 2 < liminf f \Vv R \ 2 
Jb d «->o J Bo r^o J Bb 



lim 

H->0 



l-X R 



-I \V(Uo WR ) 

b p\ b p(1-\ r ) 



< lim 

R->0 



(l-\ R ) d - 2 [ \Vw\ 2 + CUp(U) 2 f \Vw R \ 2 =( 

J Bp J Bp\B p ( 1 _x R ) JB 



\Vw\ 



Since w is arbitrary, the previous inequality implies both minimality of Uq and strong conver- 
gence U R — > Uo in W 1,2 (B p ; U{m)) as R — > and concludes the proof. □ 

The final ingredient is the following small-energy regularity result in the spirit of the fun- 
damental e— regularity theorem for harmonic maps proved in |SUlj . The result is similar to 
|Mos[ Proposition 4.1] but the stationarity condition as well as the argument of the proof 
there (the so-called "moving- frame" trick) are different. Here we modify the elementary ap- 
proach to regularity of |CWY] for harmonic maps into spheres, by rewriting the right hand 
side of (|3.14p - (|4.3p in a suitable way and applying a standard estimate for the Laplace equa- 
tion. Then, an iteration argument gives the decay of the BMO norm at small scales, whence 
continuity follows from the equivalence of Morrey-Campanato spaces and Holder spaces in a 
suitable range of parameters. 

Proposition 4.11. Let d > 3 and U G W 1,2 (T^;U(m)) be a weak solution to the equations 
(I3.14P . Then there exist e > and (3 > 0, both independent of U and ko £ T^, such that if 



sup sup R 

keB RQ (k o )0<R<Ro 



2-d 



\VU\ 2 dk < e 



B R (k) 



for some R > 0, then U £ C°^(B Ro/2 {k );U(m)). 

Before going into the proof we quote two auxiliary results. First, recall that by definition 
for an open set Slcl^a function u £ (fl) is in BMO(Q) if 



(4.21) 



where up 



sup -f \u 
D Jd 



\U\\BMO(Q) = 

u is the average of u over D and the supremum is taken over all balls Dcfi. 



D 



The first fact we need is a classical result of John and Nirenberg (see [St], Chapter 4). 
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Lemma 4.12. For any 1 < p < oo, there exists a constant C p > (which depends only on p 
and d) such that if u G BMO(£l) then 

(4.22) \\u\\bmo(Q) < 

where the supremum is taken over all balls D <Z Q, C M. d . 

The second auxiliary result is a standard regularity property for solutions to the Laplace 
equation. 

Lemma 4.13. Let d > 3 and C M d &e an open ball of radius R > 0. Lei g G (^-,2), 
s = There exist C > depending only on q such that if F G L 2 (i?^;]R d ), 5 G L 2 (B R ) 

and u G Wo ,2 (-Br) is a u>eafc solution to An = divT + i/ien 
(4-23) ||V«|| L9(Bs) < C (H^IU^Bfi) + • 

Proof of Proposition \j^Tl First note that if [/ satisfies the condition 



sup sup R 2 ~ d / |W| 2 d£; < e for some i?o > , 

k£B RQ (k o )0<R<R o JB R (k) 
then on Br q = Bn (ko), by Cauchy-Schwartz and Poincare inequality we have 

(4.24) 

\\U\\bmo(b R() ) < sup [i \V-i U\ 2 ) <C sup (R 2 i \VU\ 2 ) <Ce 1 ' 2 . 
d r cb Ro \Jd r Jd r J d r cb Ro \ Jd r J 

Now, we aim to show that, for e sufficiently small, there is a quantitative decay of the BMO 
norm of U at smaller and smaller scales. 

Up to translation we may assume ko = 0. Let a G (0, |] a fixed number to be specified 

later. For each k G Br /2 an d ^ G (0,-Ro/2], let Z?t = Dt(k) C I?r be an open ball of 

radius i, and for each k G D at = D at {k) let -R G (0,i) be such that B a ji(k) C -Do-t- Clearly, 

JB H /jfc 5 C D t C B Ro , so that if R G (R/2,R) we still have the bound -R 2_d / \VU\ 2 dk < e. 

JB R (k) 

On the other hand, given a constant matrix To G M m (C) with |To| < e -S- ^0 — 

$B R (k) ^ 1 we ma y choose .R G (R/2,R) so that 

(4.25) / \U-T \<8[ \U-T \. 

JdB R (k) JB R (k) 

Since f]aB s G W 1 / 2 ' 2 \dB f>;U(m)) there exists a harmonic extension /i G W l,2 (Bf i ; M m (C) so 
that h=C/|gs^ on d-B^ and /i G C°° in the interior. Moreover, mean value formula, Jensen's 
inequality and (|4.25p easily give 



(4.26) \Vh{k)\ p <C P R~ P £ |*7-T | P 

JB R {k) 

for any p G (l,oo) and any k G Bi^(k). 
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On the other hand, as in the two dimensional case, since U £ W 1,2 (T^', U(m)) is a weak 
solution to (|3,14p . we have 



(4.27) 



j=l 3 1 



for some 1? function / such that U 1 f £ L 2 (T^;u(m)). As in Proposition 14. 3| if we set B J 



1 T j*du au* 

2 \ u dkj dkj 
is a weak solution to 



\ U *W] ~ l?kJ U ) then B3 G L ' 2 ( T *d> u ( m )) for 3 = 1 , . . . , d and U — h€z W^ 2 (B R (k); M m (C)) 



d r)TT d fl - /f) R 

(4.28) A(U-h) = ^^ + f = ^£((U-T )B^) + f-(U-T )^ dBl 



■ ^ dk 3 
3=1 J 



1 dk 3 



3=1 3 

Since \B{k)\ < C\VU(k)\ and \g{k)\ < C\f(k)\ < C(l + |W(A;)|), for the right hand side 
of (|4.28p we have the straightforward estimates 



(4.29) \\F\\ Lq{B m < \\VU\\ LHB Ck)) \\U-T \\ ^ < \\U - T \\ * y/iR 



d-2 
2 



and 



(4-30) \\9\\Ls(B R (k)) < CR° + C\\VU\\ LS(Bnm < CR^ (R 2+d ^ V + y/iR 



Applying Lemma [4. 131 to (I4.28P and taking (14. 29ft and (I4.30P into account we obtain 
(4.31) 



q - „d-2 



B R (k) 



\V(U-h)\ q < C q eiR q ^\\U -T c " 9 



Oil 2q t 
L^(B R {k)) 



i.e. 

(4.32) 



/ \V(U-h)\ q < C q R~ q ei I -f 

JB^{k) \JE 



, 2q 

\U-To\*=i 



2-q 
2 



+ C a R 2q + C a e 1 2R qd2 -^ +q -2- d 



<B k {k) 



Now we choose p = q* = -jj^- > q and R = aR < R < R. Using Sobolev inequality, ()4.32p 



and (I4.26P we estimate 



B<r R {k) 



c_ 

~R d JB h{ k) 



\U-h{k)\?<^- d \ p-hf + ^l _\h-h(k)\* 



c_ 

R d j 



< 



p 

C ^r (V \ v ( u ~ h ^ q ) q + C ^ P su p \ Vh \ p 

R \JB n (k) J B m 
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V 

2-q 

2 



<C<j~ d [ei(i \U-T Q \^\ +C q R 2q + C q e q 2R qdl ^ +q i- d ) +Ca p i \U-T^ P 

\ \JB s (k) J I JB R (k) 

+ R 2p + ef R pd^+P§- P -d \ +C(T p I p _ To | P 

/ JB R (k) 



2-q 

\U-T \ — 

B R {k) I I JB R (k) 



Since Bn{k) C Dt, if we choose To = f BR ^\ U the John-Nirenberg inequality (|4,22p yields 
(4.33) 

(f \U- h(k)\ p ) < U-^e 1 ' 2 + a) C q \\U\\ BMO m t )+C q a- d / p (t 2p + e§ t pd ^ +p ^ p - d 

\JB aR (k) J 

On the other hand, by Holder inequality and (14.33P we get 

/ \ 1/2 / \ l/p 

(4.34) / _\U--f U\<lf _ \U-h(k)\ 2 ) <(/ _ \U-h(k)\ p ) 



Bcr R (k) JB aR {k) \JB aR (k) I \JB aR {k) 



< (a-d/P £ Vi + ^ C q \\U\\ BMO{Dt) + C q a~ d / p (t 2p + £ h pd -^ +p i- p - d 

hence, taking the supremum over B a n(k) C D a t we obtain 
(4.35) 

\\U\\ B MO { D at) < (a- d l p e'l 2 + a) C q \\U\\ BMO{Dt) + C q a~ d l p (t 2p + e f ^+*>f-^ 

Since p = p(q), s = s(q), the exponent a := pd 2 -^- + p^ — p — d — > 1 as q \ so 
we can fix q small such that a G (0,2). If we choose a G (0, |] and e > so small that 
C q {a- d / p e^ 2 + <r) < \ then 

(4.36) ||f/||BMO(D CTi ) < \\\U\\ B MO{D t ) + Ct a , Vt G (0,i?o/2] , 

where C > and a G (0, 2) are independent of U, k and t. Thus, from (|4.36p an elementary 
iteration argument on v(t) = ||t^||,BA/0(D t ) an d the John-Nirenberg inequality (|4.22p give 

(4.37) / ^W-i ^U\ 2 < Ct 2f} , Vt G (0, Rq/2] , Vfc G B Ro/2 , 

JDt(k) JDt(k) 

for some f3 = /3(a) > 0. Finally, from |Ca| we conclude that U G C®'P{Bn Q /2\lA{m)) and the 
proof is complete. □ 

Combining the previous propositions and the Liouville type theorem proved in the Appendix 
we have the main result of this section. 

Theorem 4.14. Let d = 3 and U G W 1 ' 2 (T^ l ;U(m)) be a minimizer of the localization func- 
tional ([33]). If 2 < m < 3 then U G C°(T*;W(m)). 

Proof. Fix feo G Tg and for each R = R n \ we define Un{k) = U{k$ + Rk). According 
to Proposition 14.61 such a sequence is bounded in W^ 2 {M?;lA{m)) and, up to subsequences, 
it converges weakly to a degree-zero homogeneous map U Q G W{j t (R 3 ;U(m)). According to 
Proposition I4.8[ such convergence is strong and the limiting map Uq is a degree-zero homo- 
geneous local minimizer of the Dirichlet integral in W^ 2 (M. s ;U(m)). According to Corollary 
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IA.101 when 2 < m < 3 we have U {k) = const, therefore ^- J Br , & , \VU\ 2 = J Bi \VU R J 2 ->■ 
as n — > oo, hence it can be made arbitrarily small at sufficiently small scale. Thus, in view 
of Proposition 14.111 continuity around ko follows, and U G C°(T^;U(m)) since ko was arbi- 
trary. □ 



5. Analytic regularity 



In this section we first prove analytic regularity for continuous weak solutions to the Euler 
Lagrange equations (j3. 14|) . whenever x is a real-analytic Bloch frame (i.e. Aj G C""(Tj, u(m))). 
Then, combining this property with the continuity results for the minimizers of the localization 
functional fj3 - 5f> . we prove analyticity for any minimizer of the functionals (13. 4p and (13. 5p . 

We start with the following auxiliary result. 

Proposition 5.1. Assumed > 2. Let \ be a real-analytic Bloch frame andU G W 1,2 (T d ; U(m)) 
be a weak solution to equation (l3~T4"p . If U e C°(T* d ;U(m)) then U G J^ 1 ' 4 ^; W(m)) and 
C/ G W 2 > 2 {Y* d -U{m))- 

Proof. We rewrite the system (|3.14p in the form 

(5.1) AU = J r (A(k), U, VU) = F(k, U, VU) , 

where T : (M m (C)) d x M m (C) x (M m (C)) rf — > M m (C) is a real-analytic (polynomial) map and 
the matrices A{k) = (Ai(k), . . . , Ad(k)), i.e. the entries of the Berry connection given in (|3.6|) . 
are real-analytic on the torus T^. It is easy to see that, since A is smooth and U takes values 
into U(m), by construction the function F(k,s,p) on the range of U satisfies the structural 
assumptions 

(5.2) \F(k,s,p)\+\V s F(k,s,p)\ <c (l+|p| 2 ), |V fc F(fc,s,p)|+|V p F(fc,a,p)| < ci(l+|p|) 
on T* x U(m) x (M m (C)) d . 

As a consequence of |Joj . Lemma 8.5.1 and Lemma 8.5.3, any continuous weak solution U of 
(15. ip is locally in W 2 ' 2 H W 1 ' 4 and the conclusion follows taking a finite cover of the torus. □ 

Combining the previous result with the standard regularity theory for linear elliptic equa- 
tions and the fundamental analyticity results for nonlinear elliptic systems (see e.g. |Mor| . 
Chapter VI), we obtain full regularity. The proof is standard, so we just sketch it for the 
reader's convenience. 

Proposition 5.2. Let d = 2 or d = 3 and let U G W 1,2 (Y d ]U(m)) be a weak solution to 
equation (I3.14p . If x is real-analytic and U G C°(T d ;lI(m)), then U is real- analytic. 

Proof. According to Proposition Owe know that U G W l ' A {Y* d M(m)) and U G W 2 ' 2 (T* d ; U{m)). 
If d = 2 then Sobolev embedding yields VU G IP for any p < oo. Similarly, in case d = 3 we 
have S7U G L 6 , hence (j5.2|) implies G(k,U,S/U) G I?. Thus, linear elliptic regularity for (|5.ip 
gives ?7 G VF 2 ' 3 and in turn VU £ L p for any p < oo again by Sobolev embedding (note that 
the same property holds for any d > 4, compare |Jo[ Lemma 8.5.4]). Clearly, if VU G IP for 
any p < oo the same is true for G(k, U, VU) because of (15. 2p . hence linear elliptic regularity 
for (|5.ip gives £7 G VF 2 ' 2 ' for any p < oo, which in turn yields U G C 1,a for any a G (0, 1) 
by Sobolev-Morrey embedding. Going back to (I3.14P a standard bootstrap argument in the 
Holder spaces C l ' a , I > 1, yields by induction U G C l > a => AU G C l ~^ a => f7 G C* +1 ' Q , so 
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that U G C°°(T^; U(m)). Finally, since the coefficients Aj in (|3. 14j) are analytic, by the results 
in |Morj . Chapter VI, any smooth solution is real-analytic. □ 

The main result of the section is the following. 

Theorem 5.3. Let 1 < d < 2 and m > 1, or d = 3 and 1 < m < 3, or d> 4 and m = 1. Let 
X be a real- analytic Block frame and Fmv(') an d Fmv{' \ X) th> e Junctionals defined by 
and \3. 5\) respectively. Then: 

(i) any minimizer U G W l,2 (T* d ; U(m)) of Fmv(' > x) * s real-analytic; 

(ii) any minimizer ip = {<fi, ■ ■ ■ , y m } C W of Fmv(') is a real- analytic map from M. d to 
(U { ) m . 

Proof, (i) Since any minimizer U is a weak solution to (|3.14p and \ ls real- analytic, the 
conclusion follows from Proposition 14.11 if d = 1 or m = 1, and from Corollary 14.41 Theorem 
14.141 and Proposition 15.21 in the other cases. 

(ii) Let w = Uftpcp and ip = x ■ U. Then Fmv{w) = Fmv^) = Fmv{U ; x), is a minimizer 
of the latter functional in view of (|3.7p . and the conclusion follows from part (i) above. □ 

6. Exponential localization of maximally localized Wannier functions 

The main result of the paper is the following, and provides an affirmative answer to problems 
(MVi) and (MV 2 ). 

Theorem 6.1. Let er* be a family of m Block bands for the operator U.l\) satisfying the 
gap condition \2.9(l . and let {P*(k)}k£R d be the corresponding family of spectral projectors. 
Assume d < 2 and m > 1, or d > 1 and m = 1, or d = 3 and 1 < m < 3. Then there exist 
composite Wannier functions {wi, . . . ,w m } C W which minimize the localization functional 
(|3.2p under the constraint that the corresponding quasi-Bloch functions are an orthonormal 
basis for Ran P*{k) for each k G Y* . In addition, for any system of maximally localized 
composite Wannier function w = {wi, . . . ,w m } there exists /3 > such that e^ x ^w a is in 
L 2 (M. d ) for every a G {l,...,m}, i.e. the composite Wannier function w a is exponentially 
localized. 

Conjecturally, we expect that the parameter (3 appearing in the latter claim does not depend 
on the minimizer w, and that the claim holds true for any (3 < a, where a is appearing in 
(l2~10j) . 

Proof. In view of Theorem 12.91 there exists a Bloch frame x which is real- analytic. Therefore, 
problem (MVi) is equivalent to showing that the r.h.s. of (|3.7p is attained, which is proved in 
Theorem 13.41 

Let w = {w\, . . . w m } C W be any minimizer of Fmv- Then ip a := U-QFW a defines a 
minimizer {<fi, ■ ■ ■ , (p m } C W of Fmv among the Bloch frames. With respect to the frame x, 
one has if = x • U, where U G W 1,2 (T* d ; U{m)) is given by Ub ja (k) = (Xb(k) , ip a (k)). Clearly, 
U is a minimizer of Fmv{' > x)- By Theorem l5.3l U G W 1,2 (Y d ;U{m)) is actually real-analytic. 

The function U defines a real-analytic r*-periodic function U : M. d — > U{m). By unique 
continuation, U extends to a holomorphic function U c from to GL(m, C), which is T*- 
periodic in the real part of its argument. Analogously, arguing as in the proof of Proposition 
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12, 1| x admits a holomorphic extension x C £ ^r/3 2 ^ or some P2 > 0. Therefore 99 = % • U is 
a real-analytic Bloch frame which admits a holomorphic extension <p C = X C ' U c , which is in 
'K^ o Q for /5g = min{/3i, /32, a}. Moreover, for any (3 < (3q there exists C such that 

\\<p c (k + ih)\\ 2 dk < C 

Y* 

for every h such that \hj\ < (3 for j G {1, . . . , d}. By Proposition !2.3l one has that := ip a 
satisfies 

y e 2/3|x| |w a (»| 2 dx < +00. 

□ 
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Appendix A. Harmonic maps into U(m) 



We consider, for Q' = R 3 \ {0} and U G W^(Cl'; U{m)) with m > 2, the energy functional 



r 1 3 / f)TJ* f)TT \ 

(a.i, E(c , ;Si) = y n _g tr (__j dt , nccn , 

We assume that U is a local minimizer of (jA.ip in O', i.e. that E(U;£l) < E(W; O) for 
any CC ft' and for any W G W^ft'; such that supp(C7 - W) CC ft. Clearly, if 

* G Co°(ft; u(m)) and eel, then U E (k) = U(k) expe^(fc) is an admissible variation of J7, 
hence local minimality gives 



E(U E ;n) = 0, -J 

e=0 ae 



£(f7 £ ;ft) > 0. 



£=0 

Since the tangential variation ^ can be chosen arbitrarily, the first condition in (|A.2j) easily 
implies that U is a weakly harmonic map, i.e. U is a weak solution to 



loc 



We aim to prove that, when m > 2, any local minimizer [/ G ^^(O'; U{m)) which is 



degree-zero homogeneous (a minimizing tangent map), i.e. 

U(k) = uj f for some u G ^(.S 2 ; W(m)), 

is constant. The argument we are going to use combines a stability inequality derived from 
(|A.2p (see inequality ()A.12p below) and a nontrivial quantization property for the energy of 
every harmonic map u) G C°°(S 2 ; U{m)) |Val Corollary 8]. 

Actually, in view of Lemma I A. II below, when m = 2 this constancy property is known and 
it follows from |SU2l Proposition 1], since SU{2) = S 3 (\^2). However, when the target is a 
sphere, the constants in the stability inequalities are uniformly bounded as the dimension of 
the sphere increases (see |SU2| formula (*)]), which will not be the case in the problem we 
are dealing with. Here we prove, by a different technique, the constancy property in the case 
m < 3. In our opinion, if (|A. 12[) cannot be improved, then it seems difficult to prove the 
Liouville property for any m > 4 using the so-called Bochner method as in the sphere-valued 
case (see |SU2| and |LW2| ; see also |Xin| Chapter 5] and references therein). 

As far as our specific problem is concerned, we can assume that the target is indeed the 
special unitary group SU(m) in view of the following auxiliary result. 

Lemma A.l. Let d = 3 and m > 2. Let U G W^(ft'; U(m)) be a degree-zero homogeneous 
weakly harmonic map. Then det U = a G and U G WjJ?(Sl'; SlA{m)) up to multiplication 
by a constant unitary matrix Uq. As a consequence, tr(U* djU) = for each 1 < j < 3. 

Proof. Clearly, det U G W^(fl'; U(l)) and it is de gree-zero homogeneous. According to [BZJ, 
if B is the unit ball, there exists g G W 1,2 (B) such that det U = e 19 a.e. in B. By slicing, 

1 2 

g is Wj ' on a.e. ray from the origin, so it is continuous along a.e. ray. Since e ig = det U is 
constant along the rays, we conclude that g is also constant along the rays, i.e. g is degree-zero 
homogeneous. 
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Let us set U(k) = e™ 9(fc) I and W(k) = U(k)*U(k). By construction det U = detU, so 
that W = U*U £ W 1,2 (B; SU{m)). Thus, in order to prove the lemma it is clearly enough to 
show that g (and in turn U) is constant in B, because the conclusion follows in the whole M 3 
since U is degree-zero homogeneous. Notice that if 77 £ Cq°(B), e£M and g e = g + er\ then 
(A.4) 



E&*IW;B) _ /^|V*|» + /J|> * = E(ei'-I; B) + E(W;B) . 

Since U e = e™ 9s IW is an admissible variation for U, differentiating (|A.4j) we readily see that 
the function g is weakly harmonic in B, hence g is continuous (real-analytic). Since g is also 
degree-zero homogeneous we conclude that g is constant in B as claimed. As a consequence, 
U is constant, det U is also constant (both in B and in K 3 , both functions being degree-zero 
homogeneous) and U € (^'j SU(m)) up to multiplying by a constant unitary matrix 
Uq = U*. Finally, since U*djU G su(m) we also have tr (U*djU) = for each 1 < j < 3. □ 

A.l. The stability inequality. Throughout this section we assume that U is a smooth 
harmonic map and we denote by V the variational vector field, i.e. V(k) = j^\ £ _ U e (k) = 
U e (k)fy(k) associated to the deformation U e (k) = U(k) exp £ ^( k \ We regard SU(m) C M m (C) 
as a Riemannian manifold with the metric induced by the embedding in M m (C), the latter 
being equipped with the Hilbert-Schmidt inner producO^ (A, B) = Retr(A*5). The second 
variation formula for the energy [LVVTl Chapter 1] yields 

E(U £ ;n) = Y] / \\V ej V\\ 2 - RidUie^dUiej)^) 

A JO- 



< A - 5 ' 



s=0 j 



where {ej} is an orthonormal basis of M 3 = TO', V is the pull-back via U of the Levi-Civita 
connection on TSU(m), R is the curvature (4, 0)-tensor on TSU(m), and ||j4|| = (tr(A* A)) 1 / 2 
is the Hilbert-Schmidt norm in M m (C). 

In order to obtain a convenient stability inequality, we focus on a variational vector field V 
which is obtained by orthogonal projection onto TSU(m) of a constant vector field on M m (C) 
and we average over such constant vector fields the corresponding inequality given by (|A.5j) . 
The idea is not new, it is explicitly used in |SU2j when the target is a sphere (averaging over 
the conformal vector fields) and more generally in |HW| for homogeneous manifolds and e.g. 
in |We| for general targets. Here we follow a very concrete and elementary approach when the 
target is SU(m) C M m (C) and we obtain a stability inequality with an explicit constant (see 
inequality (|A.12j) below). 

Remark A. 2. Following |HW| the same inequality (with exactly the same constant) could be 
deduced in a more abstract way, regarding SU{m) as a homogeneous (group) manifold. More 
precisely, one can regard SU{m) as minimal submanifolds in the sphere S m ~ l {y/m) through 
the standard minimal immersion in the first nontrivial eigenspace of its Laplace-Beltrami 
operator, getting a stability inequality with an explicit constant expressed in terms of the 
first nonzero eigenvalue of the Laplacian (see e.g. |Xin| . pages 137-138 and equation (5.23) 
or |HW| . pages 328-329 and Proposition 5.2). The eigenvalue as well as the dimension of the 



( 2 ) Notice that this metric differs by a constant from the metric on SU (m) induced by the Killing form on 
su(m). The difference is, for our purposes, immaterial. 
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eigenspace are known in the literature, in terms of the representation theory of the Lie algebra 
su(m), so the constant can be explicitly computed. 

To implement this idea, we consider the orthogonal projection Pjj : M m (C) —> TuSlA{m) 
defined by 

PjjU) = -u (u*(f) -<ft*u - — ti(u*<t> - 4>*u) i 

2 \ m 

Clearly, for £7 = 1, the formula above reduces to the orthogonal projection from M m (C) onto 
su(m). For a fixed (ft £ M m (C), we define the tangent vector field (ft J along U by setting 

4>T(k) = P u{k) (<p). 

For T) G Co°(n',M), we define V^^ik) := r]{k)ft{k), as a section of the pull-back bundle 
U*TSU(m), corresponding to the admissible variation uf ,r, (k) = U(k)ex.p(eri(k)U*(k)(ft~ [ (k)). 
The advantage of this choice is that the second variation of the energy, when averaged with 
respect to (ft varying in an orthonormal basis (ONB) of TM m (C) = M m (C), decouples as the 
sum of two simpler terms, as shown in the following lemma. 



Lemma A. 3. Let V^'* 1 be defined as above, and llf ,ri be the corresponding variation. Then 
E(U+";Sl) = (m 2 -l) f \Vr,\ 2 



V — 



de 2 

fliGONB 



e 
e=0 



(A.6) + [ V 2 Y1 E {||V ei T || 2 -i?(dC/( ei ),0 T ,dC/( ei ),^)} 

3 0eONB 

where the sum runs over (ft varying in an orthonormal basis of M m (C). 
Proof. By the Leibniz property, for V = V^ ,v = n(ft J one gets 

l|V ej ^|| 2 = \d jV \ 2 \\^\\ 2 + 27]^ (ft, V e ^ T )+// 2 ||V e ^ T 

We notice that 

(ft ,4> J )= Yl ^' PjJ ft =^ p U = dimSU(m) 

</>£ONB 0GONB 

where Tr denotes the trace in the algebra End(T[/M m (C) J3EJ Moreover, 

(ft , Ve 3 ft) =dj \ ' ^) I = d i (\ ^SU{m) \ = 0, 

</>eoNB y <f> / \ ' 

so one obtains 

2 El|V ej ^|| 2 = dimcSW(m)|V77| 2 + 7? 2 ^ ||V ej .0 T || 2 . 

0gONB j 0GONB j 

By substituting in (|A.5P and recalling that dimSU(m) = m 2 — 1 one obtains the claim. □ 



In contrast, all over the paper we denote by tr the trace in M m (C), i.e. the ordinary matrix trace. 
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We now exploit the specific structure of SU{m) to make the last term in ()A.6[) more explicit. 
The first lemma does not depend on the particular structure of V r t >,v , so we state it for any 
variational vector field V. Hereafter, we set Cj = and, in view of the embedding SU(m) C 
M m (C), we identify dU(ej) and djU. 

Lemma A. 4. Let V = Utp with tfi(k) £ su(m) and U(k) G SU{m). Then 
(A.7) l|V e3 T/|| 2 - R(djU, V, djU, V) = \\djip\\ 2 + tr (U*djU [ip, djijA) . 

Proof. The covariant derivative in the pull-back bundle U*TSU(m) is equal to the projection 
on U*TSU{m) of the ordinary derivative, i.e. 

V ej V(k) = P u(k) (|£(*o)- 

By an explicit computation, and taking into account that t/j* = —if), one gets 

v ej v{k) = u(k) (?L(k) + \ [w(k)d 3 u{k),m\ 



Thus, one directly computes 



(a.8) \\v ej v\y = WdrfY+ixixrdjU - a,-^*)) + - \\[u*djU, v]||- 

Since SU{m) is a Lie group with bi-invariant metric g(A, B) = Retr(j4* B) for A,B€ su(m), 
the curvature tensor R can be written in terms of Lie brackets by the Cartan formula 

(A.9) R(A, B, A,B) = l tr ([[A, B] , A]*B) = I \\ [A, B]\\ 2 . 

Thus, by left invariance, one has 

R(djU, V, djU, V) = R(U* djU, ib, U* djU, if>) = - || [U*djU, i/t] \\ 2 
which cancels exactly the last term in (|A.8j) . yielding the claim. □ 
Lemma A. 5. With the definitions above, one has 



(A. 10) £ £ {llVe^ll 2 -^^^^^)^.!^^. 

3 </>eONB j 



U\\ 2 . 



Proof. By Lemma [A. 41 applied to V(k) = 4> J (k), one immediately gets 

||V ej <A T || 2 - R(djU,ft JjUrf) = \\d^f +tr (U*djU[iP, drf]) 
where ib(k) = U(k)*(p J (k). By setting 

$ = -(U*6-<j>*U\ i.e. ib = $ tr(^) I, 

2 m 

one obtains djib = djip — — tr(djib) I. Since djtp and I are orthogonal in M m (C), 

\\d^\\ 2 = \\d^\\ 2 --\tr(d^)\ 2 . 
m 
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We first prove that 



E ^|tr(^)| 2 = l||^C/f. 



m m 
</>eONB 



Indeed, by the Parseval lemma 

|tr(9^)| 2 = ^(djU^-rmi 2 



4 

0GONB </>€ONB 



£ \\2BBtr (M-idjU))? 



4 

0SONB 

= \{4>,-id j U)\ 2 = \\d j Uf. 

0GONB 

The sum of the remaining terms vanishes, after the summation over an ON basis. Indeed, 
since [ip, djip] = [tp, djip], one has 

tT(U*djU[i/}, djip]) = \\dj^\\ 2 - -lv{(U*djU<j)*U + U*<j>U*d j U)(d j U*<l>-(l)*d j U)} . 
Therefore, 



{||^|| 2 -tr(^C7[V, drf])} 

fliGONB 

(A.ll) = jti(2U*d j U<l>*Ud j U*<l>-U*d j U<i>*U<l>*d j U + U*<j>U*d j UdjU*( l 



0SONB 

We sum over <f> in the canonical basis {$> a p, ^aft} > for a, f3 € {1, ... , m}, where 

[®a/3] ab = Saafap and = i$ aj3 . 

Every term in (jA.lip in which only (ft (resp. only <fi*) appears, provides a null contribution to 
the sum over <j), since the term containing <£ Q/ g cancels the one containing & a p. Therefore, 

{ ll^f-tr (iTdjU [if,, drf])} = ^J2 2tl ( U * d i U ®«,P Ud i U *^) 

<j>eONB a,fi 

= \ti(U*djU)\ 2 = 0, 

where we used that U takes values in SU(m). □ 

By the previous lemmas, we obtain the following conclusion. 

Proposition A. 6 (Stability inequality). Let U E C°° (fl' , SU (m)) be a local minimizer of 
IA . 1\) in SI' . Then, for every Q CC ft' and every rj € Cq°(J),M) one has 

,2 II tt— 1 jj-7-||2 ™^™2 i\ / |V7„|2 



(A.12) / rf \\ir L dU\y < m{m A - 1) / \Vr][ 

Jn Jo, 

Proof. We consider the variation Us corresponding to the variational vector field V(k) 
r](k)(j) J (k); by minimality one gets 

d 2 



de 2 



E(U^; 0) > 0. 

e=0 
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By summing over <f> in an ONB of M m (C), and taking into account Lemmas I A, 31 lA~4l and 
IA.5l one obtains 

o < £ $2 E{ut^ si) = {ra 2 - 1) f \vn\ 2 - 1 f r? £ ll^t/|| 2 , 

e=o J " ■' " j 

which by left invariance proves the claim. □ 

A. 2. Quantization of energy and constancy of tangent maps. We aim to apply the 
stability inequality (|A.12p to degree-zero homogeneous minimizing harmonic maps (the so- 
called tangent maps which appear as blow-up limits of a given minimizer when rescaled around 
a given point, as in Proposition 14, 8p . in order to show that they are constant. 

To prove this property, we restrict ourselves to the case Q' = ]R 3 \ {0} and we assume U to 
be a degree-zero homogeneous harmonic map (so that the determinant will be constant in view 
of Lemma lA.lj) . i.e. we suppose now that U(k) = uj(k/\k\) for some map uj : S 2 — > U(m). 
If U € W£ C 2 (R 3 ; U(m)) is a local minimizer of (lA.lj) . then w is a smooth harmonic map 
uj G C 00 {S 2 ;U{m)) (actually uj G C°°{S 2 ; SlA(m)) up to multiplication by a constant unitary 
matrix, in view of Lemma[AT|). Indeed, for m > 2, since U is weakly harmonic and degree-zero 
homogeneous, then uj is also weakly harmonic, i.e. it is a critical point of the energy functional 

(A.13) £(uj) = \( \\uj-\ 

2 Js 2 

defined on W 1 ' 2 (S 2 ;U(m)). Since any critical point of ()A.13p is C°°-smooth |LWll Chapter 
3], one concludes that uj is C°°-smooth, as claimed. 

Now we take degree-zero homogeneity of U into account. By localizing the inequality 
(|AT2]1 on S 2 , i.e. by taking r){k) = p(\k\)ip (k/\k\) with p G Cg°((0,oo)) and V G C°°(5 2 ) and 
optimizing in p (see |SU2l Lemma 1.3]), one obtains 



l duj\\ 2 dVol 



(A.14) f \i;\ 2 - Woj^dujlf dVol < -m(m 2 - 1) f ( |V^| 2 + -\^\ 2 ) dVol . 

Js 2 2 2 Js 2 \ 4 / 

In view of the definition ()A.13p . we set ip = 1 in the previous inequality and we derive the 
following proposition. 

Proposition A. 7. Let m > 2 and uj G C 00 (S 2 ;U(m)) an harmonic map. IfU(k) = uj 
is a local minimizer of (jA.ip . then £(uj) < —m(m 2 — 1). 

In order to proceed, we recall that a very precise description of the space of all the harmonic 
maps uj G C°°(S 2 ; U(m)) was given in |Uh| . proving the so-called factorization into unitons. 
While the latter paper exploited algebraic techniques, a slightly different factorization result 
was obtained in |Va| . based on an energy induction argument. 

Proposition A. 8 ( |Va| . Corollary 7'). Let uj : S 2 — > U(m) be a nonconstant harmonic map. 
Then there exist a natural number I > 1 and a canonical factorization 

(A.15) u = Uo(pi-pi)---(pi-pt), 81it<£(uj), 

where ujq G U(m) and each pj is the hermitian projection onto a sub-bundle of S 2 x C m 
holomorphic w.r.to the complex structure induced by the operator d + dpi + . . . dpj~i- 
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Here we regard each projection pj as a map in a complex Grassmannian m (C) , k = 
Rankpj, and each factor pj — pj- in (|A.15[) as a corresponding map into U(m), through the 
isometric embedding Gfc jm (C) U{m) (the Cartan embedding) defined by assigning to each 
subspace a unitary operator corresponding the reflection w.r.to the subspace. 

The main result in |Va| shows that each factor in (|A.15P changes the energy by an integer 
multiple of 8tt. The consequence which will be relevant to us is the following quantization 
property for a critical point of ()A.13j) . 

Proposition A. 9 ( |Va| . Corollary 8). The energy £{uj) of any harmonic map oj : S 2 — > U{m) 
is an integer multiple of8ir. In particular, if u is noncostant, then £(oj) > 8tt. 

A straightforward consequence of Propositions IA.71 and lA~9l is that £{oj) = 0, and hence oj 
is constant (and in turn U is constant) whenever ^ m(m 2 — 1) < 8ir, i.e. for m = 2. As already 
mentioned, this way we recover the regularity property for minimizing harmonic maps into S 3 
proved in [SU2l Proposition 1], since SU{2) = S 3 (V2). 

The case m = 3 requires additional care but the conclusion still holds, hence we have the 
following result. 

Corollary A.10. Let 2 < m < 3 and let U{k) = oj (|£y), U £ W^£(R 3 ; U(m)) be a local 
minimizer of (|A. 1[) . Then U is constant. 

Proof. As already recalled at the beginning of this subsection, oj € C°°(S 2 ]U(m)). Since 
the case m = 2 follows readily from Proposition IA.71 we assume m = 3, and, arguing by 
contradiction, we may assume U ^ const. By Proposition IA.71 we obtain £{oS) < 12n, hence 
S(uj) = 8ir because of Proposition IA.9I Going back to Proposition IA.8I we see that, up to a 
constant unitary matrix, we may assume that the product in (|A.15P contains only one factor. 
More precisely, since S 2 = CP 1 , Gi^C) = CP 2 and, up to isometry, Gi^C) = G2,3(<C), then 
u = p — p i , where p : CP 1 -> CP 2 is an holomorphic map (equivalently, an holomorphic line 
bundle over CP 1 ), and £{u) = 8tt. 

Now, as in |Val page 131], let Cj be the Kahler form on CP 2 normalized to be the positive 
generator in H 2 (CP 2 ;Z,). Then, in homogeneous coordinates \z\,Z2[ G CP 1 , p is a polinomial 
map of degreee one, since 8tt\ deg p\ = \f cpl p*u\ = 8tt\ci(p)\ = £{oj) = 8ir, where c\{p) 
is the first Chern number of the holomorphic bundle p. As p is nonconstant, p{z\,z?) = 
\pi(z\, Z2), P2(z±, Z2), P3(z\, z%)\ S CP 2 for three suitable (not all proportional) degree-one 
polynomials. Since the degree is one, they are linearly dependent, hence, up to a linear (resp. 
unitary) change of coordinates on CP 2 (resp. on C 3 ), we may assume p% = 0. To summarize, 
we can regard p : CP 1 — > CP 1 C CP 2 (the range of p being the set of horizontal lines in 
C 3 ) and, up to a further reflection in the third coordinate of C 3 , oj : S 2 — > U{2) C W(3) 
through the diagonal embedding (the unitary operator being now the identity on the third 
coordinate). As we are back to the case m = 2, i.e. U is a ZY(2)-valued local minimizer, we 
conclude U = const, and the contradiction concludes the proof. □ 

Combining the previous corollary with the well known regularity theory for harmonic maps 
in three dimension (see jSUlj . [Si] and |LW1] ) we have the following result which seems to be 
of independent interest. 

Theorem A.ll. Let 2 < m < 3 and £1 C M. an open set. IfU € W£ c (n;«(m)) is a local 
minimizer of (lA.lj) then U is real- analytic. 
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